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Abstract 

> 

We present several formulae for the large t asymptotics of the Riemann zeta 
function C( s )i s = a + it, < a < 1, t > 0, which are valid to all orders. A 
particular case of these results coincides with the classical results of Siegel. Using 

t— i these formulae, we derive explicit representations for the sum n ~ s f° r certain 

ranges of a and b. In addition, we present precise estimates relating this sum 
with the sum ^ c n s_1 for certain ranges of a, b, c, d. Finally, we derive certain 
novel integral representations for the basic sum characterising the leading large 

• rH t asymptotics of C( s )- 

X 
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1 Introduction 



It is well known, see for example theorem 4.11 in Titchmarsh [9], that the Riemann 
zeta function ((s), s = a + it, a, t e E, satisfies the equation 



cm = E 



l-a 



n: 5 



1 - s 



It < 



27TX 



C> 1, cr > 0, t -)■ OO. (1.1] 



Here we present a systematic investigation of asymptotic formulae for C( s )- I n this 



connection, in section we derive the following exact formula (see lemma 2.1): 



C(i 



n=l 



1 



(2tt)* 




2 s dz 
e z - 1 



< 77 < 00, -- < fa < J = 1,2, (1.2) 
where the contours of integration in the first and second integrals in the rhs of equa- 



tion (1.2) are the rays from r/exp(— in/2) to ooexp(z0x) and from ?7exp(z7r/2) to 
ooexp(i0 2 ), respectively. 

1.1 The large t asymptotics of £(s) valid to all orders 



Equation (1.2) suggests a separate analysis for the cases t < 77, t = 77 and t > 77. 



The first two cases are analysed in theorems and ^2 of section [3j Theorem JCL 
presents the large t asymptotics of ((s) valid to all orders in the case of (1 + e)t < 77. 
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For example, the formulae of theorem 31 yield the following equation for the leading 
asymptotic terms in the case when 2£ < 77: 



cm = £ 



1 



n 



n=l 



+ 



n 



1 - s V2tt 



l-s 



2ii] 



(27T) 1 " 



f JST 

-iarg(l - e 171 ) + — ReU 2 (e tv ) 







2t < 77 < 00, < a < 1, £ -> 00, 



1.3) 



where the error term is uniform for 77, a in the above ranges and the polylogarithm 
Li m (z) is defined by 



oo u 

Z 



Li m (z) = ^— , m>l. 



k=l 



Theorem 3.2 presents the large £ asymptotics of ((s) valid to all orders in the case 
of 77 = £. The formula of theorem 3.2| yields the following equation for the leading 
asymptotic terms: 



[£1 



C00=£ 



n- s - 



n=l 



l-s \2tt 



l-s 



+ 



2it~ s e u 
(2tt) 1 - s 
t~"e H 
(2n) l ~ s 



6 icx I 

- — [t + i{a - 1)] - TmLi!(e rf ) + — ReLi 2 (e rf ) - -ImLi 3 (e* 

l±l^rt_ii( 3(T _2) 

+ 24^^( fa2 - fc + 1 ) + l^( 45ff3 - 45ff2 + 4) 



1+j 

576£§ 



(36a 4 - 24a 3 - 24(T 2 + 12a + l) \ +0(r a ~ 2 ) 



< a < 1, £ — >• 00, 



;i-4) 



where the error term is uniform for a in the above range. 

The best estimate for the growth of £(s) as £ — >• 00, is not based on equation (1.1 ) 
but on the well known approximate functional equation, see for example equation 
(4.12.4) on page 79 of Titchmarsh [5]: 



n<x n<y 



x 2/ = 75 — > < (J < 1, 

Z7T 



1.5) 



where x(s) is defined by 



(2tt) s . /7rs 
X(«) = sm ( -7T 

7T V Z 



r(i- s ) 



with r(s) denoting the Gamma function. Building on Riemann's unpublished notes, 
Siegel in his classical paper [S] presented the error terms of the rhs of eq uati on ( 



4.1 



and 



.5) 



4.4 



to all orders, only in the important case of x = y = ■Jtjl'K. In theorems 
we present analogous results for any x and y valid to all orders. Theorem 4. 1| yields 
the following explicit result for the subleading term in the rhs of equation (1.5) when 
1 < 7] = 2ixy < \ft. For every e > 0, there exists an A > such that 



127T J 



n=l n=l 



Tfl 



+ e- lns T{l- s)e~^r ] a - 1 x 



e * 



1 - e -M » 



1 < i] < AU < oo, 

. I\ L 3 \ 1 r— 

0{e V + ^), U < rj < AVi < oo, 
dist(?7, 27rZ) > e, < a < 1, t ->■ oo, 



where the error term is uniform for all 77, cr in the above ranges. In corollary |4.3 
similar formulas are presented for the case 



The case of e\/t < 77 < t is analyzed in theorem 4.4, which provides a formula valid 
to all orders in terms of the function $(r, u) defined by 



$(r,«) 



^ttitx +2iriux 



-dx. 



t < 0, u e C, 



;i.6) 



0\1 



where \ 1 denotes a straight line parallel to e 37r4//4 which crosses the real axis 
between and 1. The formula of theorem 4.4| yields the following result for the 
leading asymptotic terms: For every e > 0, 



'/ -1 L Ztt j -1 

n=l n=l 



'pn — sj^ g 2 t] e v v 



a - 1 

irj 



77 

2vr 



7T7 — z?7 1 $ 



O e 2 



£5/6 



Vt<T] <t, < (7 < 1, t 



OO, 



where the error term is uniform for 77, a in the above ranges and $ and c?2$ are 
evaluated at the point 



2nt 


2t 


2nt 


r ^ 1 






-8 


9 ' 

7] z 


V 


t] 2 


L271-- 




.77. 





1.8) 



Remark 1.1 1. In the particular case of 77 = v27rt, equation (1.7) agrees with the 
analogous formula of Siegel [8]. 
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2. Although Siegel suggested already in [8] the possibility of deriving a formula 



such as (1.7) based on the function (1.6), the authors have not been able to locate 
such a formula in the literature. In the case of o = 1/2, an alternative very interesting 
rigorous asymptotic representation for £(s) involving a sum which is smoothly rather 
than sharply tru ncat ed is presented in [3]. 



3. Equation (1.7) is particularly useful in the case when t] 
a rational number. 



^ where b > is 

b 



Indeed, in this case the function $ and its derivatives evaluated 



at the point (1.8) can be computed explicitly (see equation (4.29) below). Therefore 



in this case theorem 4.4 yields an explicit asymptotic expansion of ((s) to all orders. 

4. Equation (1.3) is only useful as an asymptotic formula in the case when t/t] — > 
as t — > oo (otherwise the error term is as large as the retained terms). Therefore 
the asymptotic range where t < rj and t/r] = 0(1) is not covered by this equation. 



However, this asymptotic sector is covered by corollary 4.3 



1.2 The explicit form of certain sums 



Theorems 



3.1 



and 



3.2 



allow us to evaluate the sum Yl a n s •> f° r certain a and b, to all 



orders, see theorems 5.1 and 5.2 For example, the leading order of such sums follows 



immediately from equations (1.3|) and (|1.4|): For every e > 

1 



n 



1 - s 



2tt 



l-s 



m 

271 



l-s 



1.9) 



'1 + e)t < rji < r] 2 < oo, < a < 1, t — > oo, 



I*] 
E 



1 



n 



1 - s 



+i 



2tt 



l-s 



t 



;i.io) 



(1 + e)t < rj < oo, < a < 1, t -)■ oo, 
uniformly with respect to i]i,r]2, rj, and a. 

1.3 The explicit form of the difference of certain sums 

It is well known, see for example page 78 of Titchmarsh [9], that for large t, the 
following sums coincide to the leading order: 



.s-l 



n s and xO) ^ 

x<n<N t < ra < — 

Theorems 14.11 and 14.41 allow us to evaluate the difference of similar sums to all or- 



ders, see theorem 5.3 In particular, equation (1.7) immediately implies the following 
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equation: 



n=[— ]+l 



E ^ = x(«) E ^ 



+ e T(l-s)<e 2 77 i e" l2 ' rl L27rJ /; 



x 



$ + ^_Ji( ( )...,<]> + ( 2 
277 



2tt 



7TZ — %T\ 1 $ 



\\Tt < rjx < r] 2 < t, < a < 1, t 



— )■ 00. 



'/2 



1.4 Several representations for the basic sum 

The remarkable fact about the large t-asymptotics of the Riemann zeta function is 
that, whereas the higher order terms in the asymptotic expansion can be computed 
explicitly, the sum appearing in the leading order is transcendental. In section [6] 
we present several integral representations of this basic sum with the hope that these 
representations may be useful for the estimation of this fundamental sum (some results 
in this direction can be found in pQ). 

Let C* denote the semicircle from %r\ to it with Re z > 0: 



C77 



i{v + t) t 



71 



< 6 < 



71 



;i.i2) 



The following representations are derived in section [6] 



e 2 
(2w) s 



v s-l 



-dz 



05 



(27T) 

2 



lim / — 



t p s-l 



U 



s-1 



-du 



(27T, .,„ 

[5?] 

E 



dp + 0(t 



7? 



s-l 



<T-1\ 



< i] < t, t — y 00, 



where it is assumed that dist(r/, 27rZ) > 5, dist(t,27rZ) > 5, and the contour in the 
second equality above denotes the principal value integral with respect to the points 



{27m 


n e z + , 


r 77 j 


+ l<n< 


r t 1 


} 




1271. 




L2ttJ 





Basic Assumption: Throughout this paper it is assumed that 

77 > 0, t > 0, r}<£ 2ttZ, t £ 2ttZ. 



;i.!3) 
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2 An Exact Representation of £(s 

In this section, we assume that the branch cut for the logarithm runs along the negative 
real axis. 

Lemma 2.1 (An exact representation of Riemann's zeta function) Let ((s), 
s = a + it, a, t G R, denote the Riemann zeta function. Then, 



C(s) = X(s) J^n- 1 + 



n=l 



(27T) 



// 1 17TS 



+e" 





e 








1 - e~* 



< T] < CO 



7T , 7T 

» -o < & <o- 3 = 1,2, 



(2.1) 



where 



X ( a ) = i-L sin (— J T(l - a) 



(2.2) 



mi/i r(s), s6C, denoting the Gamma function, and where the contours of the first 
and second integrals in the rhs of (2.1) are the rays from r]exp(—in/2) to ooexp(i0i) 



and from 7] exp(iir /2) to ooexp(i02) respectively. An equivalent formula is 



ai-s) = j2n s - 1 



n=l 



(27T) 



-— + I e~ 

s 



ins J 


e z z s 1 dz 








1 - e~ z 



0<r) < oo, -| < fa < |, j = 1,2. 



(2.3) 



Proof. The Riemann function £(s) can be defined by 

v «-l 



cm 



2wr Ju a e~ 



-dz, 



s G C, 



(2.4) 



where H a denotes the so-called Hankel contour, which is the union of the following 
curves in the complex z-plane cut along the negative real axis (see figure [T]): 



H a = {ae id , -it < 6 < tt} U {\z\e 



a 



< \z\ < oo} U {|z|e l7r , a < \z\ < oo} 



< a < 2tx. 



(2.5) 



We will first derive equation (2.1) in the particular case < r\ < 2tx. In this case 
replacing r) with a we find 



CM 



X{s) 
(2tt) 



i_X I Z7TS 



e 2 / +e 2 



e z z s 1 



dz 



< a < 2tx; 
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-ix IX 

~2~ < * < 2' 



1 — e~* 
J = 1,2. 



(2.6) 



A 




>• Re 2 



Figure 1 ITie decomposition of H a into Li + L 2 + L%. 



We decompose the contour H a into the union of three different contours (see figure 
[l]), namely 

H a = L 1 UL 2 U L 3 , 



where 



7T 



< 9 < 7r| U {|^|e i7r , a < |z| < 00)}, 

e", -tt<0<-| ; , 

£<•<£}■ 

The integral along Li can be written as follows: 

— dz = ^— z ^dz = e™ s \ 6 U , (2.8) 







L 2 = {|z|e 


L 3 = < 





where the curve L\ is defined in (2.7) and the contours of the second and third integrals 
in (2.8) are the rays in the cut complex z-plane from aexp(iir/2) to ooexp(i7r) and 
from a exp(— in/ 2) to 00 respectively. 

Indeed, in the domain enclosed by L\ and by the ray from a exp(i7r/2) to 00 exp(i7r), 
i.e. in the shaded domain 1 of figure [2j we have 

z = re* 9 , a < r < 00, — < 6 < it; 



thus, 



e 



2 |= e rcose , a <r<oo, O<cos0<-1. 
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Figure 2 The domains 1,2, and 3. 



i a 



I 



-ia 



Hence, Cauchy's theorem implies the first equality of (2.8). Then, the substitution 



u 



z exp(— zV), yields the second equality in (2.6). Similarly 



v s-l 



L 2 



- 1 



-dz 



ae ~2~ 



1 - e z 



l dz 



ae^ e -u u s-l 



1 - e~ u 



-du. 



(2.9) 



where the curve L 2 is defined in (2.7) and the contours of the second and third integrals 



in (2.9) are the rays in the cut complex z-plane from ooexp(— in) to aexp(— in/2) 
and from oo to aexp(z7r/2) respectively. 



Indeed, the first equality in (2.9) is the consequence of Cauchy's theorem applied 



in the domain enclosed by Li and by the ray from oo exp(— in) to aexp(— in/2) (i.e. 



in the shaded domain 2 of figure |2J), whereas the second equality in (2.9) follows from 
the substitution u = zexp(in). 

The integral along L3 can be written as follows: 



v s-l 



-dz = — 



+ 



L 3 



' ae TT 



~ x dz — 2i sin ^ 



ns\ a 



(2-10) 



where the curve L 3 is defined in (2.7) and the contours in the first and second integrals 



in (2.10) are the rays in the cut complex z-plane from a exp(— in/2) to 00 and from 
00 to aexp(z7r/2). 
Indeed, 



v s-l 



-dz 



L 3 



L 3 1 



- l dz 



- l dz. 



(2.11) 



L 3 



The second term in the rhs of (2.11) yields the second term in the rhs of (2.10). 



Furthermore, in the domain enclosed by L3 and by the two rays (aexp(— in/2), 00) 
and (00, aexp(in /2)), i.e. in the shaded domain 3 of figure |2j we have 



re 



16 



a < r < 00, 



n n 

< 9 < -: 

2 2' 



thus 



a < r < 00, cos# > 0. 



Hence, Cauchy's theorem implies that the first term of the rhs of (2.11) equals the 



first term in the rhs of (2.10) 
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>► Rez 



Figure 3 The curves CI, CZ* and C%. 



Adding equations (2.8)-(2.10) we find equation (2.6) but with 



0. How- 



ever, Cauchy's theorem implies that the rays (a exp(— in/2), oo) and (oo, aexp(iir/2)) 
can be replaced with the rays (a exp(— in/2), oo exp(z0i)) and (ooexp(i0 2 ), a exp(i7r/2)), 



respectively, and then we find (2.6). 



In order to derive equation (2.1) with 2tt < r\ < oo, we introduce the curves C/l 



-v 



and C°, which are the following semi-circles in the complex z-plane, see figure |3 



CI 



cm 



c« 



2 + 2 ' 
~i(v + at) , (v-<x) w 

\~ „ C- 



2 2 

i(rj + a) (rj — a) 



2 6 ' 



71 71 

— < 9 < - 

2 2 

7T 7t 

— < 9 < - 

2 2 



-7T < 9 < U (| < 9 < 7T 



(2.12a) 
(2.12b) 
(2.12c) 



The second integral in the rhs of equation (2.6) can be rewritten as an integral along 



the curve C% plus an integral along the ray (^exp(z7r/2), ooexp(z0!)). Similarly, the 



third integral in the rhs of equation (2.6) can be rewritten as an integral along the 



curve —C_^ plus an integral along the ray (rj exp(i7r/2), ooexp(i0 2 )). Hence, the sum 
of the second and third integrals in the rhs of equation (2.6) yield the sum of the 



second and third integral in equation (2.1) plus the sum of the following two integrals: 
h = -e~r I —z s ~ l dz, J 2 = \ —z s ~ x dz. (2.13) 



cz 



l-e- 



v 1 — e 



Making the change of variables u = zexp(in) in the integral I\, we find 



h 



! 7T - 

2 



-u 



~ l du 



-u 



~ x du + e 



(",) 



u 



~ x du 
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17TS 

e 2 



—u^du + - 77 s ). 



(2.14) 



The sum of I2 and of the first term in the rhs of equation (2.14) can be evaluated 
using Cauchy's theorem: 



h + h = (2tt) s Y m 3 - 1 + - 



m=l 



Hence equation (2.1) follows. 



Multiplying equation (2.1) by %(1 — s) and using the identity 

X(s)x(l - s) = 1, sGC, 
together with the functional equation 

C(l-s) = x(i-s)C(s), ^C, 

we find the equation ( 2.3[ ). 



(2.15) 

(2.16) 
QED 



3 The Asymptotics of the Riemann Zeta Function 

for £ < 77 < 00 

We are mainly interested in the behavior of £(s) in the critical strip < Res < 1, 
thus we will henceforth assume that a G [0, 1]. It is possible to extend the results of 
this section to the case when a belongs to any bounded interval, but the proofs are 



more complicated (e.g. equation (3.12) involves powers of a that cannot be neglected 
ifa>l). 

Theorem 3.1 (The asymptotics to all orders for the case (1 + e)t < rj) Let((s), 
s = a + it, o, t e M., denote the Riemann zeta function. Then, 



[JLl 



CM = E 



1 



n 



V 



1 - s \2ir 



1-s 



(3.1) 



71=1 

in(l-s) 00 N-1 

:.ee- ' 

n=l j=0 



(2vr 



1 d 



'"-(l-s) 00 N-1 



+ 



+ 



e 2 



(27T 



z— it In 2 



n=l j=0 



n + ^ dz J n + - 
1 d^ j 



z 



o(^(2N + l)!!jv(- 



+ ^ dz J n+ %t 

e ^2{N+l) 



z=ir) 



Z= — lt] 



n 



-a-N 
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(1 + e)t < rj < oo, e>0, 0<(7<1, N >2, t ->• oo, 

where the error term is uniform for all rj, e, a, iV m t/ie above ranges and (2N + 1)!! is 
defined by 

(2N + 1)!! = 1-3-5 {2N - 1)(2N + 1). 

For AT = 3 equation \3.1 ) simplifies to 



fJLl 



cm = £ 



n — 



n=l 



7/ 



1 - S V27T 



1-s 



(3.2) 



+ 



2ii] 



(2tt; 



1-s 



f — ?(T 

-zarg(l - e"n + ReLUe 27 ') 

V 



+ [it 2 - 3iat -(a- l)t -ia{a+ 1)] Im Li 3 (e ir? ) J 



+ 



n3+cr 



1 + e\8 



(1 + e)t < r] < oo, e > 0, < a < 1, £ ->■ oo, 



where the error term is uniform for all r/, e, a m the above ranges and the polylogarithm 
Li m (z) is defined by 



OO L. 

Z K 



k=l 



(3-3) 



Similarly, it is straightforward for any N > 4 to derive an asymptotic formula for ((s) 
analogous to (3.2) with an error term of order 



O 



AM 



1 + e\ 2(^+1) 



Proof. Suppose first that < t < rj < oo, e > 0, < a < 1, iV > 2. All error terms 
of the form O(-) will be uniform with respect to r],e,a,N unless otherwise specified. 
The proof is based on equation (2.3), i.e. 

r_2_i 

c(i-*) = £ 



n=l 



n- 1 - - ( ± 
s \2tt 



(J^) S + G L {t,a; V ) + G u {t,a;r ] ), < rj < oo, (3.4) 



where Gl and Gu are defined by 

,ooe*i e -z z s-l dz 



G L (t,G-;r]) 



and 



Gu(t,a;r]) 



e 2 



(27T) 



177" 5 

e 2 
(2tt) s 



1 -e- 



hi 



e z z s l dz 
l-e~ z '' 



71 71 

2 Y 2' 



71 , 71 

2 < ^ < 2- 



0<r/<oo, (3.5) 



< 7] < oo. (3.6) 
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The asymptotics of Gl 

Using the expansion 



1 oo 

- J — = Y' 

1 - e~ z ^ 

n=0 



Re z > 0, 



(3.7) 



in the definition (|3.5|) of Gl, we find 

oo 

E 



i •> -«? 



^<0i<^, 0<??<oo. (3.8) 



The interchange of the integration and the summation in (3.8) is allowed because the 



sum on the rhs of (3.7) converges absolutely and uniformly in Re 2; > 8 for any 8 > 0, 



and because by assumption 77 ^ 2ttX. Integration by parts shows that the integral in 



(3.8) can be written as 



e -nz z s- 



hi 



N-l 

- l dz= ^ e ~nz+it\nz 
3=0 

■it] 



1 d \ J z' 



3 „o-\ 



n — - dz 1 n 



nz+itlnz D N dz, 



Z = —IT] 



(3.9) 



where = D^(z,n,cr,t) is short-hand notation for 

N 



D N : 



d 1 

dz n — - 



„<T-1 



(3.10) 



We claim that 



e 2 



E 



-nz+it In 2 



Avcfe = 0^(2iV-l)!!(JV + l)V 



l-N 



(3.11) 



In order to derive equation (|3.11|) we first note that we can write Dn in the form 

(3.12) 



N N 



(AO a-l 



t b (nz) N - b a c 



z 



6=0 c=0 



(wz - it) 2N 



where A^ , 6 = 0, . . . , N, c = 0, . . . , N, are integers which satisfy 

14^1 < (2JV- 1)!! = 1 -3-5 (2JV-1), b = 0,...,N, c = 0,...,N. 



Indeed, if -Da? can be written as in (3.12), then the computation 



df \_ {N) t\nzf^ 
* bc {nz-itf N 



dz \ n — 
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A (N)„-t t b (nz) N - b a c 
be z ( nz _ it y(N+i) 



[it(b -N)- ita - nz(l + b + N) + nza)] 



shows that -Dtv+i also can be written as in (3.12) 



In order to estimate the lhs of (3.11), we choose <pi = 0. Then, since t < \z\ and 
In — — I > n for all z on the contour, we can estimate 



D H = o(vs-m<r l ±±\($ 



6=0 c=0 



(nz) N a c 



z 2N (n - ^) 2N 



■0 



0( (2A^ - 1)!!|^| CT - 1 - JV (A^ + 1) 2 t^ 



n 1 



z = u — irj, u G [0, oo) 



Thus, the lhs of (3.11) is 



OO -| 

nzjt log z | ( 0A j _ i^l/' JV_L ^ 2 \^\ (T - 1 - N fa 



e- nRez \z\ a - 1 - N dz 



Oie-^Yl / \e- nz e ltlogz \(2N -1)\\(N + lj 2 \z 

n=l 

0((2iV-l)!!(iV + l) 2 f;-l / 

^ 71=1 

( °° -| /"OO 

(2iV- l)!!(iV + l)V~ 1_7V ^^vy e ~™^ M 

o((2N-mN+vy-i- N f2j^j 



rv 



OU2N - 1)\\(N + 1) 2 t] 



2a-l-N 



This proves (3.11 ). 



Equations (3.8), (3.9), and (3.11) imply that Gl satisfies 



iir s OO N — 1 

(■l.U.rr:!!) 

n=l j=0 



(27T) 



1 d V z^ 1 



n — - dz / n 

z 



(3.13) 



z=—ir\ 



+ (2iV-l)!!(A^ + l) 2 i7 



2„o--l-JV 



t < 1] < OO, t — ^ OO. 



The asymptotics of GV 

We now let e > and suppose that < (1 + e)t < rj < oo. In analogy with equations 



(3.8) and (3.9) we now find 



Gu{t,<y;rj) 



ITYS OO 

e 2 



(2*) 



i?2 



71 



71 



e- nz z s ~ Y dz, -tt<02<77, 0<r?<oo, (3.14) 
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and 



pa 
J ir\ 



e i<t>2 



e -nz z s- 



N-l 

- 1 dz= J2e~ nz+inn 

3=0 

J ITj 



1 d V z a ~ l 



n — -dz ) n — li 

z 

nz+ithlz D N dz, 



z=ir] 



(3.15) 



where is defined by (3.10). We claim that 



27V S OO 

e 2 



(2w) 



E 



""- + '"" '-n , i - o( (2 \ ... \)]].\( ' ' ' V iV ;/" V 



=1 



z D N dz = 0\^2N- 1)\\N{ 



(3.16) 



Indeed, since can be written in the form given by (3.12), Jordan's lemma implies 
that we may choose 2 — 7r/2. Then, using the estimate 



D N = o((2N- l)!!^!"- 1 ^^ 

^ b=0 c=0 



£ \ (nz) a 



n 



N 



o((2N-l)\\\zr 1 - N (N + l) 2 {n _ t)2Nj , 



z = iX, A G [rj, 00), 



we find that the lhs of (3.16) is 



( °° rioo \ 
e ^J2. \e~ nz+itXogz D N \dz\ 
n=l ^ ' 

( 00 /-ioo 
e f]T/ e-^(2iV-l)!!|^r- 1 - iV (iV+l) 
n=l ^ 

(poo 00 
(2iV - 1)!!(JV + l) 2 / \ a ~ 1 ~ N d\ V 



n 



N 



[n 



L)2N 



dz 



N 



1 \2N 



Ol (2iV-l)!!(iV + l) 2 V 



a-N 



a 



N 



1 + e\ 2 ^ 



o((2N-l)\\N 



In the third equality above we have used the following estimates: 

,N 1 



E 

n=l 



rr 



n 



1 ^2iV 



(1 

o 



1 \2JV 
1+eJ n =2 



^ e JV(lnn-21n(n- I ^)) 



^l±e\ 2JV + ^ eJV(lnn _ 21nri+ _2_ } \ 

' n=2 ' 



(3.17) 
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o 



1 + e\2W 



iV > 2, 



where the last step is a consequence of the fact that e 1 -^ < — for e > 0. This proves 



(3.16) 



We next claim that (3.16) can be modified as follows: 



17VS OO 

e 2 



(27T) 



E 



- nz+itlnz D N dz 



(3.18) 



l + e\2(iv+i) 



0( (2iV + l)!!A^^^J"" "V" V - J 



Indeed, integration by parts shows that the lhs of (3.16) equals 



e 2 



(2*)' 



z+it In ; 



-D 



n=l 



n — 



N 



+ 



17TS OO 

e 2 



Z=IT] 



E 

n=l 



i0 2 



-712+it In z 



D^ + idz. 



By (3.16), the second term is 

»((2iV + 1)!!(JV + 1) (I^^V""- 1 



01 



while, in view of (3.12) and (3.17), the first term satisfies 

N N 



ITTS OO 

e 2 



(2tt>' 



(jv),. , (T _ 1 t 6 (m^) Ar - 6 (7 c 



71=1 



7? fe=o e=0 



JV AT 



t n 



(mr/ — it) 2N 

N-b 



6=0 c=0 



^ ri b (n- i) 

n=l ' v r? y 

(00 
^-i-^(2iV-l)!!(^ + l) 2 ^ 
n=l 



102JV+1 



rr 



1_02JV+1 



1+6' 
/ I _|_ F \ 2JV+1 

0[ v a - 1 - N (2N-l)\\(N + l) 2 ' X 



This proves (3.18 ). 



Equations (3.14), (3.15), and (3.18) imply that Gjj satisfies 



17TS OO JV — 1 



G v (t,cr; rj) 



e 2 



(2x) 



EE' 

71=1 J=Q 



-nz+it In 2 



1 d V z^" 1 



+ 0^(2jV + l)!!iv( 



n — - dz I n — 

z 

1 + e\2(iv+i) 



Z=IT) 



V 



a-N-1 



(3.19) 
(1 + e)t < 77 < 00. 
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Proof of (pO) 



Substituting the expressions (|3.13|) and (|3.19|) for Gl and into (|3.4|), we find 

c(i-«) = E 



n" 1 - ' 



n=l 



S V27T 



+ 



iTTS OO N—l 

e 2 



(2 -)- EE' 

V / n =l j=o 
"» oo N—l 

e EE' 



-nz+it In z 



-nz+it In 2 



1 d Y z 



ii — dz j n 
1 d V z a ~ l 



3 -er-l 



(2 



n=l j=0 



n 



it 



+ 0^(2iV + l)!!iv(E±f j 



1 + e\2(JV+i) 



r-jV-1 



n 



it 



Z= — lt] 



Z=IT) 



(1 + e)t < r\ < oo, i — > oo. 

Replacing o by 1 — a and taking the complex conjugate of the resulting equation, we 
find ([311). 



Proof of (3.2) 



Letting iV = 3 in (|3.1|), we find 

\JL] 
>-2tt> 

cm = E 



n 



1 / rj 



n=l 



1 - s V2tt 



i7r(l — s) 

e 2 



1-5 



-irjn 



n=l 



11,1 i(q((T - 1) + a) 
n 1 + q rjn 2 (1 + q) 3 



+ 



1 -gV-l) 2 + g(-2a 2 + q + 2)-g(g + l) 
r] 2 n 3 



Z7r(l — s) 

e 2^ 



+ (2^ ( -'" r E 



71=1 



(i + g) 5 

11,1 i(q(<r - 1) - *) 
nl — q rjn 2 (1 — g) 3 



+ o(V T ~ 3 ( 



1 g 2 (a- l) 2 + g(-2 ( T 2 + a + 2) + ( r( ( 7 + l) 
?7 2 rz 3 



(i - qY 



l + e 



7177 



The above sums can be expanded as series in q with coefficients expressible in terms 
of the polylogarithms Li m (z) defined by (3.3). For example, 

k oo 



E 



o'i'fn 



^— ' nil — q) ^— ' 

n=l v ^ 7 n=l 



E n E 



fc=0 



7177 



k oo 



fc=0 



E(^) E 



D vqn 



71=1 



7? 



fe+1 



fc=0 



E L Wi(e-) 
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where the interchange of the two sums can be justified by separating the terms where 
k = and noticing that the remaining double sum is absolutely convergent. More 
precisely, 

N ■ M , x b / N N ■ M x 



lim lim f — ) = lim [ h lim 

n^oo M-^oo — ' n — ' V nil J n^-oo \ z — ' n m-s>oo z — ' n z — ' \ nn 

n=l k=0 v ' 7 \n=l n=l k=l x 

00 e i-qn °° f t \ k 

Lii(e ) + ^ ^( n7? ) 

n=l fc=l v /7 

00 00 / j. \ 

fe=l n=l x '/ 



00 00 / , \ fe jjjn 



fc=0 

Similarly, for any j > and l,m> 1, 



fc=0 v 

Letting 77 — > —n in this equation, we find 



This gives the following expression for the Riemann zeta function: 



J -I 

n=l 

_ >7i-(i-s) r 00 / j- X fc • °°/r,o\/j-\ fc 



•fc=o x ,y ' fc=o 

fc+i 



"""s('r)(-;)"-*-' 



' k=0 



2<r 2 -a-2^,/* + 4\/ <\ t+1 T . , 
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A;=0 

<7(<7 + l)^/jfe + 4\/* X * 



' fc=0 



sen© 



k=0 

k+2 



k=0 

1 + e\8 N 



+ o(V ct - 3 (- 

By including only the first few terms in the sums over k, we find an expansion for ((s) 
in powers of K For example, including only the terms of order larger than 0(^J+^) we 
have 



rj 



cw=E»--t^( 



i_ /^_y- s 

1 - S V27T, 



n=l 

+ ^. 2 (e-0 - ^Li3(e-') 
+ ^We-) - ^Li3(e-) 

+ (^{ L '' (£, " ) % L ' 2(e "' )+ (^) 2Li3(e ' ,) 
- ^Li 2 (e") - ^Li 3 (e 



rj rj 
+ i ^Li,(e».) - ^fiiLi3(e-) 

^(^( i3+ (^)' 
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After simplification we find 



rjLl 
L271-J 



sv ; ^ l-s V2tt 



n=l 

+ 7^~\I~7 {^Im Lii (e iri ) + -Re Li^e 4 ") + i ( -) ImLi 3 (e^) 

--ReLi 2 (e l ")-^ImLi 3 (e^) 

(a - 1)£ . • ia(a + l) . , . , 
- ^^ImLi 3 (e ^) — — -ImLi 3 (e"J 



+ 



-,3+a 



t 3 + 



If 



The identity Lix (e ir ?) = - log(l - e iri ) now yields (3.2). 



QED 



Theorem 3.2 (The asymptotic expansion to all orders for the case t = 77) Let 

C{s), s = a + it, a, t G M, denote the Riemann zeta function. Then, 



cw=E 



1 / t 



?? 



n=l 



1 - s V 2tt 



»7t(1-s) OO iV— 1 



+ 



e 2 



n=l j=0 



l-s 



-712— it In 2 



(3.20) 



(2vr 

tTr(l-a) oo JV-l 



1 d 



v 7 n=2 j=0 



■nz—it In j 



n + - dz J n + 



1 d 



n + - dz J n + 



z=it 



z=—it 



+ 



t 

2^ 



l-s 2AT 



e « ^ c fc (l _g£X^y + o ( (2iV + 1)!! jVe^* 1 )*-*"" 



fc=0 



t 2 



< a < 1, N >2, t ->■ 00, 
where the error term is uniform for all a, N in the above ranges and the coefficients 



Cfc(cr) are given by equation (3.30) below. The first few of the Ck's are given by 

c o(o-) = — 
c x {a) 
c 2 (a) 



l + i 



12 



(6a 2 - 6a + 1) 



135 



(-45a 3 + 90a 2 - 45a + 4) 



(3.21) 
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C 4 (cr) 
C 5 0) 



i - 1 

432 

i 

5670 
1+j 
194400 



(36a 4 - 120a 3 + 120a 2 - 36a + l) , 
(189a 5 - 945a 4 + 1575a 3 - 987a 2 + 168a + 8) , 

(1080a 6 - 7560a 5 + 18900a 4 - 20160a 3 + 8190a 2 - 450a - 139) 



For N = 3 equation (3.20) simplifies to 



cw=£ 



n=l 



1 - s V 2vr 



(3.22) 



|^{|^[-* 2 - it(a - 1) + (a - I) 2 ] - ilmLijCe*) + ^HeU a (e«) 



it + a + a 



2 + 3a, 



3i, 



+ ~ ' t2 ' " zImLi 3 (e rt ) + ^-^ReLi 4 (e if ) + ^ImLi 5 (e lt ; 

t- s e u f 1 + i r- 1 , 

t — r— <^ Vvrt- -i(3<7 - 2) 

(2vr) 1 - s \ 2 3 V ; 

+ ^^v 7 ^ (6a 2 - 6a + 1) + — — (45a 3 - 45a 2 + 4) 



24v^ 
1 + i 

576tl 

i 



135* 

(36a 4 - 24a 3 - 24a 2 + 12a + l) 



+ (189a 5 - 315a 3 + 42a 2 + 84a - 8) 

+ 1 - % -. v/^(l080a 6 + 1080a 5 - 2700a 4 - 1440a 3 
103680*5 v 

+ 1710a 2 + 270a -139) J +0(t _<r ~ 3 ), < a < 1, t 



— > 00, 



where the error term is uniform for all a in the above range and the polylogarithm 
hi m (z), m > 1, is defined by (3. Sty . Similarly, it is straightforward for any N > 4 to 
derive an asymptotic formula for ((s) analogous to (3.22) with an error term of order 

o(t-°- N ). 



Proof. Setting 77 = t in (3.4), we find 



f-M 

1 / t \ s 

C(l-s) = ^n s - 1 --f— J +G L (t,a;t) + Gu{t,a;t), < t < 00. (3.23) 

71=1 ^ ' 



Equation (3.13) is valid also when rj — t and gives the asymptotics of Gl(£, a; t). On 



the other hand, using the expansion (3.7), we write 

G u (t,a;t) = G { l }\t,a;t) + G^(t,a;t), 
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where 



G^(t,a;t) 



e 2 



e z z s dz, 



|<02<|, 0<t<oo (3.24) 



and 



_vws oo />ooe il ^ > 2 

G<g\t,a;t) = ^-J2 e-^z^dz, ~ < 2 < |, < t < oo. (3.25) 



The asymptotics of can be found using integration by parts, whereas the asymp- 
totics of G$ will be computed by considering the critical point at z = it. 



(2) 

The asymptotics of G 



u 



Repeating the steps that led to (3.19) but with the sum over n only going from 2 to 



oo and with t = r], we find the following analog of ( |3.19 ): 

1 d V z a ~ l 



ITTS OO N—l 



G ( u'(t,a;t) 



(27T) 



EE 

n=2 j=0 



-nz+it In z 



n — - dz i n — 



z=it 



+ 0^(2N + l)\\Ne 2( - N+1 h°- N - 1 



(3.26) 



The asymptotics of G[/ 

Letting in the definition (3.24) of Gy, 02 — and 



we find 



Letting 



z = it + tp, p G [0, oo) 



l-./..-"' / e -i(/t»-ilog(l-ip))| 




1 - ipf-'dp. 



v = p — i log(l — «p), 



(3.27) 



we find 



G 



(i) 



2tt 



«- / e-" (1 -'" W) V 



P(W) 



where p(v) is defined by inverting (3.27) and 7 denotes the image of the contour [0, 00) 



under (3.27). In order to ascertain that the value of p{v) is well-defined, we consider 
in detail the map defined by 



: C \ [— i, —ioo) — > C, 
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p I—)- v = p — i log(l — ip), 

where [— z, —too) is a branch cut and the principal branch is chosen for the logarithm, 
i.e. 

v = p — zln |1 — ip\ + arg(l — ip), arg(l — ip) G (— 7T, 7r). (3.28) 
The function satisfies 

0(p) = O iff p = 0; 

0'(p) = -^ = O iff p = 0; 
z + p 

v = ^p) = - l -^ + 0(p 3 ), p^O. 

Moreover, we claim that <fi maps the first quadrant of the complex p-plane bijectively 
onto the region delimited by the contour 7 = 0([O, 00)) and the positive imaginary axis 
in the complex f-plane, see figure |4j Indeed, it is clear that <p maps [0, ioo) bijectively 
onto [0, zoo). Writing p = p\ + z'p 2 , we have 

<9Re</> \p\ 2 + P2 <9Re0 pi 



dpi p? + (p 2 + l) 2 ' dp 2 p? + (p 2 + l) 2 ' 

This shows that each level curve of the function Re in the first quadrant intersects 
the positive p-axis in a unique point and that the vector 

(- Pl ,|p| 2 + p 2 ) (3.29) 

is tangent to the level curve passing through p. Since 

dlmcp —pi <9Im0 |p| 2 + p 2 

~dpT = p 2 + (p 2 + i) 2 ' "^T = p? + (p 2 + i) 2 ' 



we infer that the derivative of Im0 in the direction of (3.29) equals |p| 2 , showing that 
Im0 is strictly increasing along each level curve. This proves the claim. 

Thus, if v belongs to the region delimited by the contour 7 = 0([O, 00)) and the 
positive imaginary axis, we may define p(v) as the unique inverse image of v under <p 
which belongs to the first quadrant. Using analyticity to deform the contour 7 to the 
positive real axis, we obtain 

" W J p(v) 
We claim that there exist coefficients {cfc(o")}^ _1 such that 

(i-ip(v)y 1 " V " J 



p(v) 



(S^c k {a)v k ' 2 + r N (yfi)\ v>0, (3.30) 

\ h — n / 



*fc=0 

23 




where r^^y/v) satisfies 



-^ N ^ v) -dv = 0{^^). (3.31) 



Indeed, since <fi has a double zero at p = 0, there exists a neighborhood V of p = in 
the complex p-plane and a Riemann surface £ which is a two-sheeted cover of V with 
local parameter A := y/v, such that the map p t— >■ is bijective and holomorphic 
V — > S. The function is analytic away from [— i, —zoo) and 0'(p) ^ for all p ^ 0, 
thus the Riemann surface E and the map \/v t— > p can be analytically extended as 
long as 0" 1 (f ) PI [— i, —ioo) = 0. The function </> maps the segments [—i + 0, — ioo + 0) 
and [—i — 0, —ioo — 0) onto the lines Ret> = — it and Rev = n respectively, hence 
yfv i — y p is well-defined and analytic for all \v\ < ir. Thus, the map 

(1 - ip) a _ e CTlog ( 1 -^ 
j '. A i — y A — A 

P P 



is holomorphic from the open disk {0 < |A| < y/n} C S to C. This implies that (3.30) 
holds with 

f(*)(0) 

c fc = ^p, k = 0,...,N-l, 
and with a function r/v(A) satisfying 

MA) = £ ^p^ fc , < |A| < \Fk. 

k=N 



In order to prove (3.31), it is sufficient to show that there exists a constant C inde- 



pendent of N such that 

\tn(Vv)\ < Cvi, v>0, N>2. (3.32) 
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Cauchy's estimates imply that the numbers 



y (fc) (Q) 



are uniformly bounded for all 
k > 0. Thus, equations (13.281) and (13.301) imply that such a C exists for v > 1. On 



fe! 



the other hand, an application of Cauchy's estimates in a disk of radius |A| = r < y/n 
also yields 



k=N V ' V y 



N 



|A| < r, 



where M r = maxui =r |/(A)|. Choosing any r e (1, y/n) we see that such a C also exists 



for < v < 1. This proves (3.32) and hence also (3.30) and (3.31). 



Equation (3.30) together with the identity 

r(a + i) 



e' tv v a dv 



t a+l 



a > — 1, 



imply the following asymptotic expansion of G^: 



JV-l 



fe=0 



fc+i 

t 2 



+ 



r 



■N+V 



N+l 
t 2 



(3.33) 



Proof of (3.20) 



into (3.23) yields 



Equations (3.13), (3.26), and (3.33) give the asymptotics of Gl and Gy. Substitution 



n=l x ' 



™ oo TV— 1 



(27T) 



e 2 
(2^ 



+ 



+, 27' 6 



EE 

n=l j=0 
oo N-1 

n=2 j=0 
27V 

it " 



1 d \ 3 z a ~ x 



-nz+it In 2 



n — - dz ) n — 

z 



1 d V -z 17 " 1 



n - - dz J n — 



z=—it 



z=it 



^^^^ + oY(2iV + l)\\Ne^ N+1 h^ N -A 
k=o t 2 ^ ' 



where we have replaced N by 2iV+l in ( |3.33[ ) and used that r(JV+l) = N\ < (2JV+1)!! 
to eliminate one of the error terms. Replacing a by 1 — o and taking the complex 



conjugate of the resulting equation, we find (3.20). 



Proof of (3.22) 



Letting N = 3 in (3.20) and using the expressions in (3.21) for {c^jg, we find that the 
term on the rhs of (3.20) involving the c^'s yields the term involving the second curly 
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bracket on the rhs of (3.22). On the other hand, using that 

z l ~ a (naz + i(a - l)t) 



d z- a 



n + - dz n + - 



and 



1 d 

n + - dz 



{nz + it) 3 

l -° (n 2 a(a + l)z 2 + in (2a 2 - o - 2) tz - (a - l)H 2 ) 



n + - (nz + it) 5 

long but straightforward computations show that the two terms involving the double 
sums on the rhs of (3.20) yield the term involving the first curly bracket on the rhs of 
([3^221. This proves l3~22l). 



4 The Asymptotics of the Riemann Zeta Function 

for < r] < t 

In this s ectio n, we consider the asymptotics of ((s) as t — > 00 with < n < t. 

and its corollary treat the cases e < n < \ft and 2ix\ft < n < —t under 



4.1 



Theorem 

the assumption that dist(r/, 27rZ) > e for some e > 0. The case when rj is of th e same 
order as y/i and the case when dist(7], 27rZ) — > are covered by theorem 
corollary. 



4.4 



and its 



Theorem 4.1 (The asymptotics to all orders for the case 1 < 77 < v^) For ev- 
ery e > 0, there exists a constant A > such that 



W 1 1 



(4.1) 



n=l 



71=1 



Tfl 



rjv-i, 

2 J ,(2fc), 



27T2 



s, P -([^]+i)%f (s-i) s -i f y" ^ '(ty jkfzy 



2\ *+4 



2 / 1 



+ e- i7r T(l-s)e-^r/ 



(T-l 



X 



o((¥) T *i- 



O iVe ^ + 



e < n < U < 00, 1 < iV < 



At 

„3 ) 



N+l 



, *3 < 7/ < y/t < OO, 1 < iV < 



At 
„a ? 



dist(?7, 2ttZ) > e, < a < 1, t ->■ 00, 

where the error terms are uniform for all rj, er, iV within the above ranges, and the 
function <p(z) is defined by 



(p(z) 
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For N = 3 equation (4-D simplifies to 

ft j m x 



(4.2) 



n=l 



n=l 



e CTs r(l - g) e -([|]+i)^ e f s -i e f 



+ 



2vri 
1 



27T 7] 



(1 - e"* 



1 - e~ ir > y/i 
(1 - c -*i)2((_^[i]_i( ff -!)+*)« + ff -1) 



77^ 

2(i- e -*»)(_ t? [i]_<( a -i) + t) 



+ 1 



iyfirrf 



+ e"^T(l - s)e"^r/ 



(7-1 



0(f), e<r/<t5<cx), 3r/ 3 < At, 

0(e~^ + £), U <7] < y/i <oo, 3r] 2 <At, 
dist(?7, 2irZ) > e, < a < 1, t -»■ 00, 



where the error term is uniform for all rj, a in the above ranges. 



Remark 4.2 It will be shown below (see equations (4.20) and (4.21)) that 



o 



m 



0(1), n > 0, (not uniformly in n). 



Thus, the term in (4.1) which involves ip( n '(0) is of order 



rf 1 



0[e—T(l-s)e . ^ 
o(e~ i7r T(l -s)e-% rj 



rr- ( 2jl 



n + 1 
2 \ 2 



e < 77 < £3 , 



(not uniformly in n) 



It follows that, as expected, the first term left out of the sum is smaller than the error 
term. 



Proof of theorem 4-1 We start with the expression for ((s) given on page 82 of 
Titchmarsh [9]: 



\JL] 

l 2 77 J 



c(*) = E- + x( s )E — + 

t — ' n s z — 4 n s 



1 e"^T(l - s) 



n=l 



71=1 



where 



2ttz 



s— 1 — - N 

- 1 



^Ij, 0<T]<t, 



(4.3) 



3=1 
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»4ni 



• 2ixi 







> Reu> 



C 3 v 



• — 2717 



• —A-ki 



• — 67ri 



C 4 



Figure 5 The critical point if] and the integration contours {Cj}\ in the complex 
w-plane. 



and the contours Ci, C2, C3, C4 are the straight lines joining 00, cr\ + ir](l + c), —cr] + 
ir)(l — c), — cr\ — (2[j-]+l)7ri, 00, where c is an absolute constant, < c < 1/2, see figure 
5 We henceforth let c = 2~ 3 / 2 so that the total length of the contour C2 is 77. The 
branch cut for the logarithm now runs along the positive real axis. The contributions 
to from I\, I3, and I4 are exponentially small by the arguments on pages 82-83 
of [9] (in these arguments, it is assumed that x = t/rj < y = i]/2tt, however this 
assumption is not used for the derivations of these estimates and therefore we do not 
need to make it here.) 

Let e > be given and suppose that 2e < r\ < < a < 1, m — [t/rj], 
dist(?7, 27rZ) > 2e, and N > 1. All error terms of the form O(-) will be uniform with 
respect to r],a,N (but not with respect to e). In order to analyze the integral I2, we 
write 




where 



(s-l)log(l+^)--^2-A^ 2 + (--[-])2 



Defining (p(z) by 



z = 



(s-l)l0g(l+A)-^2 



it 



(4.4) 
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we have 



(p(z) 



e z _ e -iT) 



We write for z in a neighborhood of zero, 



N-l 



n=0 



n=0 



Then 



s N (z) 



z N r {p(w)dw 
2ni J r w N (w — z) 



(4.5) 



where T is a contour which encircles and z but none of the poles of ip. We split the 
contour C 2 as follows: 

C 2 = C*2 U C*2, 

where C\ denotes the segment of C 2 of length 2e which consists of the points within 
a distance e from irj. We write 

h = r 2 + r 2 , 

where 



/ 9 e = e" 



-(m+lW^\s-l 



e 2 ' 2 '" 1 ip(w — ii])dw 



and 



^ = e -(m+l) iJ7 ^ S -l 



e 2 ^ <^(w — ir])dw. 



cr, 



We claim that 



I? =0[ e 2 -^=e 24 " 



(4.6) 



Indeed, the change of variables w = irj + gives 



i£ = e-Mj 



+ / e~^ A \(p(Xe^)\dX 



There exists a constant A > such thatQ 

.4 



g\ 77 L 77 J/ 



< 



z = Ae * , e < |A| < 4- 



(4.7) 



1 Here and below ^4 denotes a strictly positive constant which may change within a computation. 
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Moreover, the definition of 4>(w) implies 



\og(j)(w) — (s — 1) log ( 1 



(a-l)log 1 



W\ t it 2 

— W -w 

ir/ J 7] 2r/ 2 

w\ it ^ (-i) k - 1 

k + 2 



17] J T] z 



k=l 



Hence, for \w\ < |r/ we have 



, , . . . . . 8 t l2 1 \w\ 5 
Relog0H < |cr _ 1| log - + — \w\ - — 



and so 



|0(Ae *4 )| < e 



| CT -l|ln| + -^A 3 



£ <|A|<f 



Equations (4.7) and (4.9) imply that 



|^(Aex)|=o(e^ |A|3 ), e<|A|<|. 



This yields 



I r 2 =0 e-M? 



o-l / 2 * A 2 + ^4r|A| 3 



Splitting the integrand as 



and noting that 



we find 



___L\2 t \2 t \2 

g 12j7^ — g 24t; 2 X g 24?7^ 



2 e "^ A2 dA = of— r . 



2J = 01 e-^rf- L e ^ 



Vi 



This proves (4.6) 



We next claim that 



AT-1 



n=0 
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+ 



i< < 4* 



In order to establish (4.10), we need to estimate the error term 



In view of the assumption dist(rp 2nZ) > 2e, we have 

0(1), 



e z _ e -i V 



z\ < e. 



(4.11) 



Let us first consider the case 1 < N < 4f. Let \z\ < ||e (so that < 2e) and let T 



in (4.5) be a circle with center w = and radius pjy, where 

21 



20 



Shrinking e > if necessary, we may assume that 2e < ^rj. Then, since <p(w) is 



analytic in the disk \w\ < 2e, equations ( |4.5[ ), ( |4.8[ ), and ( |4.11[ ) yield 

sat(z) = O 



2rtp N N _ N 5 -^n 



Pn - \z\ 



z\ iy p- N iy e^\ =0[\z\ N p- N N e 



, 40 
W < Yl e- 



■■■■■ 



The function p ^e^ 3 " has the minimum C^f^ )^ 3 for 



value if 



2iV?7 3 



« 5t 



r/; Pat can have this 



1/3 



21, , (2N\ 

— \z\ < r? < 2e. 

20 1 1 ~ V 5t ' 



The assumption 1 < N < 4f implies that for A sufficiently small, we have (^) 1 / 3 ^ < 
2e. Hence, letting p N = (^) 1 / 3 ^, we find 



s N (z) = O \\z 



N 



jv ■ 

5et N 3 



2A^ S 



iV < 



At , , 20/2NV 40 



'•'^irj (4 - 12) 



For |z| < |je we can also take p^ = %\z\, which yields 



.s. v (-) = O ( (|2)*e*£ ( ^'* l)3 ] = O (e^^ lzr ) , | :| < r. 



Using (|4.12|) and (|4.13|), we estimate 

t(w— irj) 2 



-(m+l)iri^y-l 



s " Sn(w — irj)dw 
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(4.13) 



(4.14) 



(T-l„ — g- 



Oirf^e 



A(f)3r, 



e 2 " 2 A 



v-\ V f 



JV 

5et 3 



o 



t ,2 , 14 t \2 



\ — 

5et \ 3 n=i / 7] 
2Nrf) 2 2 [l 



■* N 2 



>r vi , 



Of T) a ~^e 



A 2 



r 



2 

iV + 1 



t \2 

e ss^ (iA 



1 < Af < 



At 



where we have used the following estimate to find the last equality: 



L 



_ * \2 

e dX = O e 116t > 



A(f)3r? 



- * A 2 

e 116 " 2 c?A 



01 e-^M " 



A(f)3r, 



Using the expansion 




x — > oo, 



and the inequality || < 3, it follows that the right-hand side of (4.14) is 



O Tf~ L e 



3iV\ 6 7] 

~) VtJ' 



(4.15) 



At 



This proves (|4.10|) in the case when 1 < A" < ^3 



We now consider the case when tz < r] < \[t. Using (4.11) and (4.8) in the 



representation (4.5) with T a circle with center w = and radius 2 1//3 e, we find that 

(4.16) 



s N (z) = O 



\z\ N e*£ 

p 3v 3 



z < e. 



1 

Now £3 < 1] implies that e^ 3 " = 0(1), so we can estimate 
e -{m+i)iv( ir] y-i f e^ {w - iT,)2 s N (w-ir))d 



w 



O if~ L e 



Tit N-l 



O rf- l e~^T 



^ \ 2 



t 



r 



AT + 1 



JV+l 
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This completes the proof of (4.10). 



We next claim that, up to a small error term, the contour C\ in the integral in 



(4.10) can be replaced by the infinite line C' 2 , where C 2 denotes the infinite straight 
line of which C2 is a part. More precisely, we claim that there exists an A > such 
that 



-(m+l)i7j 



- (w—ir)) 2 



N-l 



b n (w — ir]) n dw 



(4.17) 



n=0 

TTt At 

01 e 2 " 2 



2e < r] < ts, 1 < iV < ; 
oC^-ig-f-^jy^ U<r]<y/t, 1<N< 



At 

„2 ■ 



In order to prove (4.17), we note that the coefficient of b n on the lhs of (|4.17|) is 

0\rf- x e-^ I r^P K \"<]\ }. 



(4.18) 



We write the integrand as 



— * A 2 — - A 2 

e x e ^ A n , 



and the second factor is steadily decreasing for A > an d so it decreases through- 



out the interval of integration provided that n < N < The term in (|4.18|) is then 



(4.19) 



We now assume that iV < At/rf. Choosing \z\ < ff(^) 3 ^, equation (4.12) yields 

21 

bet \ 3 



b n = (S n (z) - S n+1 (z))z " = O 



2nrf 



N < n > 1. 



(4.20) 



Multiplying the rhs of (4.19) by 6 n and summing from to AT — 1, we find that the 
total error is 



O I rf- x er^e 



4q- 



JL 



N-l 



n=l 



5et 
2nr] 3 



Now the function (^3)3 increases steadily up to n = -4-, so that if n < A-i where 



A < 1/e, it is of order 



Hence, if e 3 < and iV < A\, with A small enough that ^Aln-r < ^, the total 



5e 



3e A 4 



error is 



7Tt At 

Ole 2 ^ 
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This proves (4.17) when N < A\. In the case when U<r]< y/i, we instead use 



(4.16) to find 



b n = {Sn{z) - S n+1 (z))z " = 0[e 



< 31T 5 



Then, since ta < r], the total error is 



_TTt t 2 5ti° 

O I rf~ l e 2 e ^N 



TTt At 

0(r] a - 1 e 2 



This completes the proof of (4.17). 



We finally analyze the sum 



n=0 



iV-1 „ 
. e -(^ (i77 )-l^ 6ne f(n + l) / 



n > 0. 



(4.21) 



U < 7] < Vt. 



e~^ x2 \ n d\ 



n=0 
N-l 



_ e -(m+l)i V( i v y-l fe„ e f ("+!)2V(1 + (-1)' 



?1=0 



n + 1 



k=0 ' 



Together with equations (4.3), (4.6), (4.10), and (4.17), this yields (4.1). QED 



Corollary 4.3 (The asymptotics to all orders for the case 2n\Tt <r] <^t) For 
every e > 0, there exists a constant A > such that 



•I -I L Z7T J -I 

CM = £-+*(.) 



n=l 



n=l 



— s 



)r(, 



e ^ (S) ([ J. ] + 1)i 22rt <™ 

2717 



e 4 

V J 



r JV-l 
L 2 



fc=0 



(2*)J 



(4.22) 



X 



O (f ) 



AT 27Tt 



2ntl <r]<^t<oo, 1<N <^f, 



( 2 N+l \ 

Ne'^i + (^f^) 2 J , 2nVi < v < 27rt§ < oo, 1 < N < 



Arj 2 
t ' 
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distl 



2nt 



2ttZ > e, < a < 1, t ->■ oo, 



where the error terms are uniform for all r], a, N within the above ranges, and the 
function ij){z) is defined by 



e - s i og (i+|g) + ^M£] 
e z — e v 



Proof. We replace a by 1 — a in (4.1 ) and take the complex conjugate of both sides. 



We then multiply the resulting equation by x( s ) an d use the identities 



C(s) = C(s), r(s) = r( a ), x(s)=xW 

X(s)x(l - s) = 1, x(s)C(l-s) = C(s). 



(4.23) 



This yields the following equation 
# 1 [&l 1 

c(*) = x(*)£^ + £^ 

71 

+ x0 



n=l n=l 

e iT(i-«)r(s N 



(4.24) 



r JV-l 
I- 2 



X 



E^ 



fc=0 



2iri 

2k 



g 117 ' 2rj 2 L »7 J 



2 = 



(2fc)! 



(-0' 



2rf 



2\ 



r k 



x 



e < r] < U < oo, 1 < iV < 



4* 

„3 5 



OiNe + 



Nrf 



2 \ 2 



U < T] < y/t < OO, 1 < N < 



At 



dist(r/, 27rZ) > e, < a < 1, £ ->■ oo, 



Replacing r? by ^, we find (4.22) 



QED 



Theorem 4.4 (The asymptotics to all orders for the case ey/t < rj < t) For ev- 
ery e > 0, there exists a constant A > such that 



(4.25) 



71=1 



n=l 



+ e T(l-s)<e 2 77 i e" l2,rJ 2^^ /; S N (s,T]) 
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e\/i <r) <t, < a < 1, 1 < N < At, t oo, 

where the error term is uniform for all 77, a, N in the above ranges and the function 
Sn(s,ti) is defined by 



2V-1 n 



n=0 k=0 



2tt 



TXl — 17] 



(4.26) 



2nt 


2t 


2nt 


r r/ " 




-t- 




r] z 


V 


f] 2 


L2ttJ 




■V- 





The coefficients a n are determined by the recurrence formula 



it 

ii](n + l)a n+1 = (cr - n - l)a n ^a n - 2 , 

7] z 



n = 0,l,2,..., 



(4.27) 



together with the initial conditions a_2 = a_i = and oq = 1, and the function $(r, u) 
is defined by 



$(r, u) 



irirx 2 +2niux 



-dx, t < 0, u G C, 



(4.28) 



o\i 



with the contour \ 1 denoting a straight line parallel to e 3?r4//4 which crosses the real 
axis between and 1. 

In the particular case when p,q > are integers, we have 



]_ ^ ^qfD—TTiqp—2-iriqu 

x {D-l)"e 

ln=0 



(4.29) 



>/ — /Tin 2 2— 2ninu _ ( 1) ^ 



a„ — iriqp— 2iriqu P ^ 



-e 4 2 e p 

n=0 



For N = 3 equation (4-25) simplifies to 



■ '/ -J ^7T J -| 

c(^)=E- + ^ s )E — 



(4.30) 



n=l 



71=1 



+ e T(l-s)<e 2 77 s J- e 77 L 27r J L 27r j /; 



X 



* + ~~ ( <9 2 $ + ( 2 



27] 

fr-2)(<r-l) 
2r/ 2 



2tt 



7TZ — Z7? 1 $ 



^$ + 2(2 



r/ 



2tt 
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7TZ - Z77 <9 2 $ + 2 



2tt 



7ti — it] ) $ 



„, 7rt77° 

+ e" — 



iVi < n < t, < a < 1, t 



— > oo, 



where the error term is uniform for all rj, ain the above ranges and $ and its partial 



derivatives are evaluated at the point (1.8). 



Remark 4.5 1. If 77 = y^Y wnere b = p/q>Ois& rational number, then equa- 
tions (4.25)-(4.29) provide an explicit asymptotic expansion of ((s) to all orders. The 
particular case b = 1 is the famous case analyzed by Siegel in jS] using ideas from Rie- 
mann's unpublished notes. We note that, building on the work of [5], an elementary 
(but formal) derivation of the Riemann-Siegel formula is presented in [2]. 



2. The results of theorem |4.4| are convenient for studying the higher-order asymp- 
totics of C(s) in the case when 77 = constant x However, if the order of 77 is strictly 
smaller or larger than the order of the results of theorem 4.4 are less convenient, 



because in this case the asymptotics of ((s) involves the asymptotics of the sums in 



(4.29) as p and/or q tend to infinity. In this case, the alternative representation of 



the asymptotics of theorem 41, which avoids the appearance of the above sums, is 
usually more convenient. 



3. The error term in (4.25) is uniform with respect to 77, a, N. More precisely, this 
means that equation (4.25) is equivalent to the following statement: For every e > 0, 



there exist constants A > 0, C > 0, K > such that the inequality 



'/ -1 z7r j -1 



71=1 



T(l 



n=l 



s)e 2 77 i e" l2 ' rJ ^ l27V> " S N (s,r]) 



< Ce- tTS T(l - s)e 



El 

2 



s N 



TJ 



holds for all t > K and all a, N, rj such that < a < 1, 1 < N < At, t\ft < 77 < t. 

4. The main difference between the proofs of theorems 4.1 and 4.4| can be explained 
as follows. Consider the representation (4.3) of ({s). If 77 is of strictly smaller order 



than /j 1 / 2 , the main contribution to the asymptotics of £(s) comes from the part of the 
contour C2 that lies within a distance e of the critical point 777; the remaining part of 
C*2 gives a contribution which is suppressed by a factor of the form e ^ cf. Eq. 



(4.6). Thus, the proof of theorem 4.1 relies on a direct study of the integral near the 



critical point 777 using the method of steepest descent. On the other hand, in the case 
when 77 ~ the asymptotic expansion of £(s) depends on the integral along all of 
the contour G%. The contour C2 is a line segment of total length 77 centered on the 



critical point 777. In the proof of theorem 4.4, we handle this integral by relating it to 
the function $(r, u) defined in (4.28). 
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5. The coefficients a n satisfy the estimate 



a n = O 



rf 1 



(not uniformly in n); 



indeed, assuming this estimate up to n, we find 

r ti 



a n+l — O 



n+l 



o 



rn-2] 

— 2 



o 



77 



n+l 



(not uniformly in n). 



Proof of theorem 4-4 Let £ > 0, < er < 1, and m = [t/77]. As in the proof of theorem 



4.1, £(s) is given by (4.3) where the contributions from Ii,Is, and J4 are exponentially 



small, and the error terms of the form O(-) are uniform with respect to the variables 
1], a, N as £ —7- 00. It remains to analyze the integral J2. In the neighborhood of w = irj 
we have 



l)log 



w 



17] 



w — %r\ 1 f w — %r\ 
' 2 



it] 



IT] 



+ 



it 



(w - irj) + -^(w - ir]) 2 + 
2r] z 



Hence we write 



(s-1) log - 



{w-irj), 



where <j>{z) is defined by (4.4). Define {a n }^ by 

00 

m = J2 

The identity 



a n z . 



n=0 



dz 



s — 1 £ «£ 



implies 



(irj + 2) ^ na n z" 1 



If] + 2 7] 7] 



1 - (277 + z) 



n=l 





2£z\ 


\»7 


^2 y 



n=0 



Hence the coefficients a n are determined in succession by the recurrence formula (4.27) 
supplemented with the conditions a_2 = a_i = and Oq = 1. The first few coefficients 
are given by 



a 
a 3 



a 1 



a -I 
if] 



a 2 



(<r-2)(a-l) 
2r] 2 



-2t + i(a - 3)(a - 2){a - 1) 
6?7 3 
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(a - 4)(<7 - 3)(<r - 2)(<r - 1) + 2i(4a - 7)t 
Representing <p(z) in the form 

N-l 



n=0 

we find 

Mz) - 2™ y r N < i> 

where T is a contour contained in the disk of radius rj centered at the origin which 
encircles the points and z once. 

Let 1 2 1 < f 77 (so that ^\z\ < |tj) and let T be a circle with center w = and radius 
Pat, where 

21, . 3 

— \z\ < On < -77. 

20 1 1 ~ H ~ 5 ' 



Then, by (4.8), 



27TPJV ,,V-V^4\ nf\~\N--N 



r7vW = * ^ J = v 1 * J ' N < 7^ 

The function p~ N e^r r has the minimum ( ^^s ) N ^ 3 for p = (^) 1 / 3 r/; p^ can have this 
value if 

21, (2N\* 3 
20 N - 77 "5* 



Hence, 



JV \ 

5et \"5"\ Ar 27 , 20/2JVV 



For |z| < 1 77 we can also take p^ = foM' which yields 

rjv(2) = 0^^)^e^ ( ^ kl)3 J=0(e^ W2 ), |z| < |. (4.32) 
We write I2 in the form 

J 2 = II + J 2 R , (4.33) 
where /.f an d denote the integrals 

„ -(w—iri)+ J! ^(w—i'q)' 2 —mw N—l 

i%= / far 1 - — ^— J2< w - ir i) ndw ( 4 - 34 ) 
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and 



rR 
1 2 



(iri) 



s-l 



C-2 



e w - 1 



tn{w — irj)dw. 



(4.35) 



We claim that 



(4.36) 



Indeed, suppose first that \e w — 1| > A on O2. Then we have for w G O2, \w — ir)\ < e, 
the estimate 

( -— mWi 

e 



(-—m)w 

e *> 



1 



< 



.4' 



whereas for w G C 2 , |w — 217 1 > e, we have the estimates 

(i-TO-l)Re-w 



e (^-m)to 

e w - 1 



01 



e Kr > 



x.e-Rc^ 1, Rew > 0, 



^ 1 e Rc 1/1 _ 1 I J 



(4.37) 



Rew < 0. 



It follows from these estimates that the factor 6 ^ w _ 1 in the integrand of /<f is of 0(1) 

on O2. We now make the change of variables w = ir/ + Ae 1 ? and split the integral 
into two integrals; in the first integral |A| < A(y)5r/ 5 while in the second integral 



^4(^)37/ < |A| < 77/2. In the first integral, we use the estimate (4.31) of r^, whereas 



in the second integral, we use the estimate (|4.32|). This yields 

d\ + 



rR 



01 rf-^e- 1 * 



',A-'x.V ( 7)11 



e ^ X 



* \2 1 14 t i2 

e 2 ^ 2 +29 ^ rfA 



Amur, 



Steps almost identical to those leading from (4.14) to (4.15) now show that 



rR 



0( ^7 e 2 I — 



y/ij' 



The case where the contour goes near a pole gives a similar result. Indeed, if the 
contour passes near the pole at w = 2wni, say within a distance 1/2 of it, we take it 
around an arc of the circle \w — irj\ = 1. Letting w = irj + e 20 and using (4.31), the 
total contribution to 1^ from this arc is given by 



01 r/^'e" 
0(7/ 
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1 -Hi / 
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0[ if^e 
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0| " e '{-) Vt 



where we used the assumption that 77 > e\ft. This proves (4.36) 
: 1 2 ■ We claim that there 

^(w—iri)+^-i z (w—iri) 2 —rnw N—l 



We next consider 1%. We claim that there exists a constant A > such that 



(it]) 



s-l 



J2 a n(w -ir]) n dw + 0(e-^- At ), (4.38) 



n=0 



where C2 denotes the infinite straight line of which C2 is a part. Indeed, in view of 



(4.37), if we replace C2 by C' 2 , the integral multiplying a n in the expression for 1% 



changes by 



O I rf- v e 



(4.39) 



We can write the integrand as 



— * A 2 — t A 2 

e 3^ x e ^ A r 



and the second factor is steadily decreasing for A > y ^77, and so it decreases through- 
out the interval of integration provided that n < N < At, with A sufficiently small. 



The term in (4.39) is then 



0\rfe 2 e 4 " 



77 \ n 



e ^ dA = O 77 e 2 e 4j j 



(4.40) 



Also, by (4.31), choosing z such that |z| < |r(i|:) 3 ?7, we find 



a n = {r n (z) - r n+ i(z))2;" 



2n7] 3 

21 

bet \ 3 



27777 s 



ln+1 



77 > 1. 



5et 



2(77 + 1)t/ e 



■ 

3 



(4.41) 



Multiplying (4.40) by a n and summing from to N — 1, we find that the total error is 



n=l 
N-l 



bet \ 3 



271,77" 



O rf- l e 
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Now (t/n)3 increases steadily up to n = t/e, and so if n < At, where A < 1/e, it is 

0(e^ tAlogl/A ). 

Hence if N < At, with A sufficiently small, the total error is 

0(e-^- At ). 



This proves (4.38 ) 



Finally, we analyze the sum 



N-l * 

(ivr 1 E a W - 

n=0 ^ C 2 



(w—ir))+ -^tj (w—ir]) 2 —mw 



e w - 1 



(w — ir]) n dw. 



(4.42) 



The integral in (4.42) may be expressed as 



1 



where L is a line in the direction argu> = 7r/4, passing between and 2iri. This is n\ 
times the coefficient of £ n in the following expression: 



e i( 10+ 2[i]rt^)+^(i +2[^]«-^)»-[i]uH^(ti)+2[^] W i-t^ ^ 



2^2 L2ttJ 



e™ - 1 

2-n-t r 17 1 r 
rj ^T^vrJ I77J 



(4.43) 
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^(2[A ]7ri _^) 



where the function $(r, u) is defined by ( |4.28[ ), i.e 

„ir«To; 2 +27n;MX 



$(r, m) 



o\i 



1 

2tu 



dw, t < 0, «gC. 



-i 



We rewrite the expression on the rhs of (4.43) as 



2ri 2 v l 2ttJ 



2vrt 
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1=0 



k=0 



Hence the expression in (4.42) is given by 



r, ■/■ \s-l 1 2[^-}ni-it+^(2[^]m-iri) 2 „ , \ 



(4.44) 
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where Sn(s,T]) is defined in (4.26). Equations (4.36), (4.38), and (4.44) show that I2 
satisfies 

N 



I 2 = 27ri(iri) s - 1 e^ 1 ^ 



S N (s } 7])+0[ e 



t J Vi 



+0(e 



Equation (4.25) follows by substituting this expression into (4.3). 



In order to prove (|4.29|) we suppose that u G C and r < 0. We claim that $ 

(4.45) 



satisfies the two recursion relations 

37T2 

$(r, u) = $(r, u + 1) - ^Le-? (u+ ^ 2 

\/lrl 

and 



$(r, it) = I — e 



7viT—27riu 



$(t, it - r). 



Repetitive use of these two equations yields the following identities: 

37ri JV— 1 

$(r, it) = $(r, u + N) — ^= e~^ {u+n+ ^\ iV > 0, 

and 



n=0 



(4.46) 



(4.47) 



JV-l 



$( r , M) = ^ (_ 1 )n e Wr- 2 ™« + (_ 1 ^ e ^V 2 r-2^« $ ^ ) ^ — iVT ), iV > 0. 

(4.48) 



ra=0 



If r = —p/q, we apply (4.47) with N = p and (4.48) with N = q to find 



3iri p— 1 

e 4 



$(r, u) = $(r,u + p) l=^2 e 



(4.49) 



n=0 



and 



$(r, «) = ^(_i)- e ™ 2 -2™« + (_ 1 )9 e ^ 2 T-2^ (J) ( r; u + ^ ( 4 5Q ) 



n=0 



Eliminating $(r, m +p) from these two equations, we find (4.29) 



It remains to prove (4.45) and (4.46). In order to prove (4.45), we note that 



$(r, u + 1) - $(t, u) 



p 2m(u+l)x p 2iriux 

e™™ 2 : dx 



0\1 

r , 

0\1 
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T7iTX 2 +2-n-i(u+^)x^ 



Jo\i 

e-^W [ e^ 2 dx. 
Jo\i 



The identities 



e mrx dx = / e~ my 



-niy 2 d V e 



4 



o\i </o\i vl r l Vl r 



imply equation (4.45). 



The identity (4.46) is established as follows: 

r ^nirx 2 +2nxiux 

$(r,«) = / — ; -dx 

J0\1 e e 

^kitx 2 +2niux 

; —dx 

1\0 e ~~ e 

^kit{x— l) 2 +2niu(x— 1) 
1 + / t -r t ^-fix 



0\1 

l-e^" 2 ™ / - — ; — dx 

l-e™ T - 2 ™ u <$>(T,u-T) 



^KiTX 2 -\-2ixi(u—T)x 



0\1 e e 



where the second equality follows from the fact that the residue of the integrand at 
x = 0is^. ' QED 



Theorem is valid for e\ft < rj < t. As a corollary, we can find an analogous 
result valid for 2it < rj < 

Corollary 4.6 (The asymptotics to all orders for the case 2n < rj < ^) For ev- 
ery e > 0, there exists a constant A > such that 

Z / n S ^ / n L S 

n=l n=l 



+ x(s)e^T(s) le^l^- )V^ +u+ ^-^S N (l-S, —) 



2nt J \ V J 

+ O ( e 2 



t J Vi 
\/t 

2tt < 7] < — , < a < 1, 1 < N < At, t oo, 
e 

where the error terra is uniform for all r), a, N in the above ranges and Sn is given by 

KM. 
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Proof. We replace a by 1 — a in (4.25) and take the complex conjugate of both sides. 



We then multiply the resulting equation by x{ s ) an d use the identities (4.23). This 
yields the following equation 



C(s) = X (s) 



1 



Z—e n l s 

n=l 



Z7T J -| 



(4.52) 



n=l 



+ e 



* /3iV\ 6 77 



eVi<r] <t, < a < 1, 1 < N < At, t -> 00. 



Replacing 77 by — , we find (4.51) 



QED 



Remark 4.7 The results of theorem AA are convenient for studying the higher-order 
asymptotics of £(s) in the case when 77 = constant x Vr. However, if the orde r of 
77 is strictly smaller or larger than the order of \/t, the results of theorem 4.4 are 



less convenient, because in this case the asymptotics of ((s) involves the asymptotics 



of the sums in (4.29) as p and/or q tend to infinity. In this case, the alternative 



representation of the asymptotics of theorem 4.1, which avoids the appearance of the 
above sums, is usually more convenient. 



5 Consequences of the Asymptotic Formulae 



Using theorems |3.1| and |3.2| we can compute several interesting sums. 



Theorem 5.1 Define the polylogarithmLi m (z) by (3.3). The following relation holds: 
1 



E 



n 



1 - s 



2n 



l-s 



2n 



l-s 



wr(l-s) 00 N-l 
n=l j=0 



(2vr 



1 d\ 3 z~ a 



iir(l-s) oo N-l 



+ 



e 2 



n=l j=0 



nz—it In z 



(27T) 

+ o({2N + 1)\\n( 



n+ - dz 



1 d 



n 



n + - dz J n + 



irn 



z=ir\2 
-irji 



Z = -VT]1 



1 + 6X^+1) 

71 



(5.1) 



V v e 

1 + e)t < 771 < 772 < 00, e > 0, < a < 1, iV > 2, t -»■ 00, 
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where the error term is uniform for all 771, 972, e, a, N in the above ranges. 
For N = 3 equation (5.1) simplifies to 

1 



rl?2i 

£ 



7? 



1 -s 

2ir)l 



2tt 



l-s 



2tt 



l-s 



(5.2) 



+ 



t — 7(T 

-iargfl - e lVl ) + ReU 2 (e tm ) 

m 



+ \ [it 2 - 3io-t - (a - l)t - 10(0- + l)]lmLi 3 (e* m ) 1 
Vi J 



2z7/ 2 



(27T) 1 " 



f — 1(7 

-iarg(l - e^ 2 ) + ReListe"' 2 ) 

^2 



+ \ [it 2 - 3iat - {a - l)t - ia{a + l)]lmLi 3 (e"' 2 ) 1 
V2 J 



+ 



t 6 + 



1 + e 



(1 + e)t < 771 < r] 2 < 00, e > 0, < a < 1, t ->• 00, 
where the error term is uniform for all rji, 7/2, e, a m i/ie above ranges. 



Proof. Replacing in equation (3.1) 77 with 7/1 and subtracting the resulting equation 
from the equation obtained from (3.1) by replacing 7/ with 77 2 , we find (5.1). Equation 
(5.2) follows in a similar way from (3.2). QED 



Theorem 5.2 Define the polylogarithm Li m (z) by (3.3). The following relation holds: 
1 



[JL1 



1 - s 



2vr 



l-s 



t 

2^ 



l-s 



ijr(l-a) oo JV_1 

1 e 2 \^ \^ -JV.-/7 h) : 
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-nz—it In 2 
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dz J n + 



Z=IT] 



Z=—IT) 



z=—it 



t 

2v7 



l-s 2N 
Jt 
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fe+l 
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(5.3) 
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/ (2N + l)\\Ne 2 ( N+1 ) (2N + 1)\\N{^)^ N+ ^ ' 
O ' t a+N ' fjV+N 

(1 + e)t < r] < oo, e>0, < a < 1, N > 2, t ->■ oo, 



where the error term is uniform for all rj, e, a, N in the above ranges and the coefficients 
Cfc(cr) are defined in (3.30). 

For N = 3 equation (5.3) simplifies to 



L271-J 

E 



1 



n 



]+i 



1-s 
+ 



1 



1-s- 



2^ 



t 3 + 



+ 



(5.4) 



t 3 + 

(1 + e)t < 1] < oo, e > 0, < a < 1, t ->■ oo, 

where the error term is uniform for all rj, e, cr m the above ranges and {■ • ■ } denotes 
an explicit expression involving a, t, and the poll/logarithms Li m (e lt ), m — 1, ... ,5. 



Proof. Equation (5.3) follows by subtracting equation (J3.20 
ilarly, equation (5.4) follows by subtracting equation 



from equation (3.1 ). Sim- 



3.22) from equation (3.2). 



The following equation is given on page 78 of Titchmarsh [9]: 

E i~x( s ) E zo=i- 



QED 



(5.5) 



x<n<N 



t <-}7< f 



Theorem 4.4 yields a precise version of the relation (5.5 ) with an explicit error estimate. 



Theorem 5.3 For every e > 0, there exists a constant A > such that 



E ^ = x(«) E 



(5.6) 



n=[-M+l 



+ e T(l-s)<e 2 ^ s l e v^i s^v i 2w j » S N (s,r}) 



3N\ 6 



</2 



r?=r?i 



ev^ < ?7i < 772 < t, < a < 1, 1 < N < At, t ->■ oo, 

where the error term is uniform for all r]i,r]2,&,N in the above ranges and the function 
Sn(s,t}) is defined by (4- 26). 
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For N = 3 equation (5.6) simplifies to 

[ — ] f22-l 
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+ e T(l - s)< e 2 77 e" 2,r ^ 27r 
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7ri — in J <9 2 $ 
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+ 0[ e * - 
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(5.7) 



7ri — in j $ 

eVr < 771 < n 2 < t, < a < 1, £ — )■ 00, 
where the error term is uniform for all ni,n 2 ,a in the above ranges, $ is defined by 



(4- 28), and $ and its partial derivatives are evaluated at the point (1.8) 



Proof. Replacing in equation (4.25) n with 771 and subtracting the resulting equation 



from the equation obtained from (4.25) by replacing n with n 2 , we find (5.6). Equation 
(5.7) follows in a similar way from (4.30). QED 



Remark 5.4 The relation (5.5) is a particular case of (5.6) (let x = t/n 2 and iV 

t/vi)- 



6 Several Representations for the Basic Sum 

The goal of this section is to present integral representations of the basic sum nS_1 
for certain values of a and b. 

Let the contour C%, t < n, denote the semicircle from it to if] with Rez > 0. 



Splitting the contour of the second integral in the rhs of equation (1.2 ) into the contour 
C* plus the ray from it to ooexp(i0 2 ), we find 



[JLl 



C(l -s)= y // 

n=l 



s-1 



ITY S 

e 2 



-dz 



^ s(2ix) s (2-k) s J c v e z - I 
+ G L (t,a;r ] ) + G u (t,a;t), < a < 1, 0<t<r], 



(6.1) 



where Gl and Gjj are defined in ( |3.5 ) and (3.6) respectively. The term Gu(t,o~;t) was 
computed to all orders as t — > 00 in theorem 3.2, see equations (3.26) and (3.33). The 
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term G^t, a; rj) with t < r\ was computed to all orders as t — » oo in (3.13). Thus, by 



comparing equation (|6.1|) with the representation obtained in theorem 3.2, it follows 

in terms of the integral 



[57/271 



n=[t/2w]+l 



n 



that we can express the fundamental sum ^2 
appearing in the rhs of (6.1), where the error is computed to all orders. For brevity 
of presentation, we state this result only to leading order. 

Lemma 6.1 (An integral representation of the basic sum) Let , t < 11, de- 
note the semicircle from z = it to z = irj with Kez > 0. For every e > 0, 



\JL] 
L2tt> 
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n 
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tTTS 

e 2 
(2tt) s 
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s(2n) s (2tt) 
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(27T) 



log(l - e u ) + O 



1 t 

+ 



t 2-a v 2-a J ' 

l + e)t<?7<oo, < a < 1, f^oo, 



•where the error term is uniform for all rj, a in the above ranges. 



Proof. Letting a — > 1 — a in (3.22) and taking the complex conjugate of the resulting 
equation, we find 



c(i-«)= xy- 1 



n=l 
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S V27T 
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(27T) 



(27T) 



— Z / — z 

y Tit H zcr 
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-e* 1 + 2TmLi 1 (e i ')) 
1 + i 0r (6<7 2 - 6<7 + 1) 



24 



V 7 * 



+ o(t 



<7-2\ 



The lemma follows by subtracting this equation from (6.1) and using the following 
equations which follow from (3.13), (3.33), and (3.26) respectively: 



G L (t,a;ri) = - 



17] 



G ( i\t,a;t) 
G { u\t,a;t) 



(2vr 

t s ~ 
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V 7rt H icr — 

2 3 24 



1 + e)t < 77, (not uniformly in e), 
1 + i y/n (6(T 2 - 6(7 + 1) ' 



+ o(t< 7 - 2 ), 



(2tt) s 



log(l - e"**) + 0(t 



CT-2> 



QED 



It is possible to derive an alternative integral representation of the related sum 



v [t/27rj 
2^\n/2-n 



[r,/27r] + l 



n 



s-l 
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Lemma 6.2 (An alternative integral representation of the basic sum) For ev- 
ery 5 > 0, 



n 



,s-l 



(27T) 



* p*' 1 



(6.3) 



< r] < t, dist(77, 2ttZ) > 5, dist(t, 2vrZ) > 5, < ct < 1, t -> oo, 

where the error term is uniform for all 77, cr m £/ie above ranges and the contour in the 
integral denotes the principal value integral with respect to the points 



^27Tn 


n e z + , 


r 77 j 


+ l<n< 


r t 1 


} 






1271. 




L2tt. 





Proof. Let denote the semicircle from it to £77 with Re z < 0, defined in equation 
2.12c) with 77 and a replaced by t and 77 respectively, see figure [6j Let R^(s) and 



Lj'(s) denote the following: 



and 



Rl(s) 



e 2 
f2^V 



e 2 



v s-l 



e z - 1 



s e C 



(6.4) 



(2 



7T 



.5-1 



e 2 - 1 



-dz, 



seC. 



(6.5) 



Cauchy's theorem applied in the interior of the disk whose boundary is ClUC^, yields 



s-1 



seC. 



The Plemelj formulas imply 



(27T) 



dp. 



Using (6.7) to replace ill in (6.6), we find 

2 



s-1 



dp + 2Ljf. 



(6.6) 



(6.7) 



(6.8) 



Replacing the contour in (6.5) by the union of following three segments: 
{ite w I < 6 < tt/4}, {ipe^ | 77 < p < £}, {z7/e ie | < 9 < vr/4}, 
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► Re 2; 



Figure 6 The contours C l and . 



we obtain 
IP — 1 



9(s-l) t s 



LJO 



-Ate" 



-ite ie d6 



1 ef^y- 1 

v e^ e ™ - 1 



TV l 

e 4 dp 



f pi0(s— l)„s-l 

—^A—iVe ie d9 



■ (27T) 



+ J 2 + J 3 ]. 



The assumption dist(i, 27rZ) > 5 implies that 



\Ji\ < 



f e -et f a-l t a-l 



A 



tde 



A 



:i- e -x 



0(t 



Similar computations show that J 2 is exponentially small and, in view of the assump- 
tion dist(r/,2vrZ) > 5, that J 3 is O^t' 1 ). Thus, 



and the lemma follows from (6.8) 



(6.9) 
QED 



Remark 6.3 Lemma 6.2 implies that the leading behavior of ((s) is characterized by 
the following integral: 



1 p s - 1 



dp. 



Remark 6.4 The estimate (6.9) is a consequence of the fact that the integral ap- 



pearing in the definition of LP t (s) does not possess any stationary points. Indeed, 
using 

= -Ve m2 , Rez<0, 

e z - 1 ^ 



m=0 
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it follows that 



17TS 

e 2 



oo pt 



c, 



imp+it In p o— 1 



p"- L dp. 



m=0 J V 

Candidates for stationary points occur at p* = —t/m and the inequality r\ < p* < t 
implies the non-existence of any stationary points. 

The next lemma gives an alternative representation for the fundamental integral 



i?*(s) defined in (6.4). 



Lemma 6.5 Let C l denote the semicircle from ir) to it with Kez > defined by 



equation (1.12). For every 5 > 0, 

"' u s ~ x du „ 
+ 



l 



(6.10) 



0<r] <t, dist(?7, 2ttZ) > 5, dist(t, 2ttZ) > 5, < a < 1, < e < 1, t^oo, 
where the error term is uniform for all rj, a, e in the above ranges. 



Proof. Let e G (0, 1). We deform the contour C* on the lhs of (6.10) so that it consists 
of the following three pieces: 

{i7]e' id | < 6 < e}, {iue~ ie \r]<u<t}, {ite- w | < 9 < e}. 

This yields 



l TVS 
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z s dz 
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u s du 



e- i6s dd 
1 ~ lT] J < ">' " - I 



+ it s 
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o e 



ite~ 



The assumption that dist(^, 27rZ) > 5 implies that there exists an A such that \e lT]e 
1| > A for all 9 G [0, e], e G (0, 1), and r? > 0. Thus, 



-ids 



et 

<-j, 0e[0 } e] 



and so 



1 



-de 



< 



ft _ l 



At 



Similarly, because of the assumption dist(t, 27rZ) > 5, 

3 ei - 1 



e ite~ 10 _ I 



dd 



< 



At 



This proves the lemma. 



QED 
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Lemma 6.5 implies that 



„, . . 1 , f* u s ~ 1 du 
RKs) = -j— - lim 



(2tt) s e^o / e iue ~ lt - 1 



(6.11) 



which in view of (6.6) and (6.9) implies the following alternative representation for 
the basic sum: 



1 /"* ?/ s -Mi/ Mi 

-Vlim/ ,V d " = V n- x + 0(t 

L 27T J 



In fact, in the remaining part of this section we will present arguments which suggest^] 
that on the critical line the assumption et — > in (6.11 ) can be relaxed to e 2 t — > 0, so 
that the following equation is valid: 



t 1 /•* m s Mm 
= t — r~ nm / - — u ) 

v ' e 2 t-+0 v 



(6.12) 



where it is assumed that dist(r/, 27rZ) > S, rj — 0(1), dist(t, 2nZ) > 5, and a = 1/2. 
Indeed, consider the integral 



, „ _ j. , . ins i (3 

(27r)'J2*00 = e-Tr / _ 



e 2 z s dz 



Instead of the contour CL we consider the finite ray from the point it] to the point z*, 
see figure [7j where 



z* = i v + (t - r/)e i( f - £) = [ v + (t- r])e- ie }e^ , < e < 



7T 



We claim that 



y* 1 6 



(6.13) 



Indeed, the assumption dist(i, 27rZ) > 5 implies that there exists an A > indepen- 
dent of r], e such that |1 — e~ z \ > A on the contour from z* to it. Thus, using the 
parametrization z = ite~ id , < 9 < e, in the lhs of (6.13), we can prove (6.13) as 
follows: 



_ 17TS 

e t 



. 1-e" 



/7~ 



1-e 



()( f' T / e^dfl } = 0(7" o 



-lie 



< t° 



a te a 



6 e -t£ne e et de 



A 



2 The rigorous justification of some of these arguments remains open. 
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Figure 7 The contour of integration of the integral in the rhs of equation (6.14) 



where we have used that 



\9 -sin 01 < 6 s 



< 9 < 1. 



It follows from (6.13) that 



l TVS 

e 2 



Ct. 



e~ z z s ~ 1 dz _ 1 t LJL , f e~ z z s 1 

= e 2 lim / cfe. 

e->0 / if 1 — e~ z 



-Z gS—l, 

1 - e- 



(6.14) 



We next employ the following parametrization which maps z* to p = 0: 



z = r/e 2 + £ ( 1 — — — p ) < ' - 



// : / (1 - | - pj ( 



e 2 



< p < 1 



We replace in equation (6.14) e z with the expression 



g -2 _ e -ite-"+ipte- 4e -ir?(l-e- iE ) _ e it(p-l)e~ iE +0(»;e) 



e ->■ 



and we also replace z in the expression z s of equation (6.14) with 



z = te^Me'^ ( 1 - — e 
M 



(6.15) 
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where (M, /z) are defined via the equation 



1 + I - - 1 J e- ie = -Me - *", M(t, e) > 0, fi(t, e) > 0. (6.16) 
J) ) T) 



Then, the rhs of equation (6.14) is given by 



e-S.0 



1-? „it(p-l)e 



t s lim e-^Af 8 -^-^' -15 / ' 6 1N . ( 1 - — e'^ 

I — e it(p-l)e- lt 



M 



dp, (6.17a) 



where 



ij)(t,e) =e-fji. 
The functions M and \x are given by 



e ->■ 



(6.17b) 



(6.18a) 



and 



Mt,e) = e-f + 0(e 2 ), e 0. 



(6.18b) 



Indeed, consider the triangle D\ formed by the intersection of the following finite 
rays, see figure |8j 



(0,1), l,l+(3-l)e— J, (0,-Me- 



t 



J] J J \ V 
The cosine rule for the triangle D\ implies the following identities: 



-M\ =\+(t-\\ -2 (--I I <-os(7r-r) 

V J \V J \v 
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1 + 



rj 



t 



1+ ( --1 

V 



1 cose 



1 (1 - cose). 



Thus, 



or 



-M 



M 



^ 2 (i-l)(l-cose 



1-2 



7] ?7 Z 



'1 — cos e), 



t t 2 

which yields (6.18a). The sine rule for the triangle D\ implies the identity 

sin fj, sin(7r — e) 



i-1 

v 



-M 
v 



Hence, 



sin u = — — — sine, 



which, using the expression (6.18a) for M, yields (6.18b). 

Equation (6.18b) implies that ip > 0, thus, we can define (W,w) via the equation 



1 - +\e 



P_ -%i> 



1 + 



M M 
The functions W and w are given by 



We iw , W(p,t,e)>0, w(p,t,e)>0. (6.19) 



W{p,t,e) = l-p + 



e 2 7] 



O 



e 2 r] 2 



e -»■ 



and 



P 

t 1- p 



1 + 



e 2 r] 



+ 



e 2 r] 2 



e -> 0. 



(6.20a) 



(6.20b) 



Indeed, consider the triangle D 2 formed by the intersection of the following finite rays, 
see figure |9] 

(0,1), (l, 1 + -£«.<(«-*)) , (0,We™). 



M 

The cosine rule for the triangle D 2 implies the following identities: 

-,2 



w2 = 1 + w- 2 ii^' 
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1 



Figure 9 The triangle 



Hence, 



W 



M 



1 + 



2^- 



(1-*) 



r(l — COS "?/>). 



Thus, using the definition (6.17b) of ip, together with equation (6.18b), we find 



e 2 rj 2 



M ' 21-fjK t 2 



o 1 JL 
W=l P ' M 



+ 0(e 4 



e -> 0. 



Substituting in this equation the expression (6.18a) for M, we find equation (6.20a). 



The sine rule for the triangle D 2 implies the identity 

sin w sin ip 



p/M 



W 



which, using the expressions for ip (equations (6.17b) and (6.18b)), for M (equation 



(6.18a)) and for W (equation (6.20a)), yields equation (6.20b) 



Using the identity 



1 _ JLp-W 
M 



s-l 



as well as the expressions for M, fi, W and w obtained earlier, it is possible to compute 
all the expressions appearing in the rhs of ( |6.17a[ ): 



Jit+a-l)lnM 



M s 1 = e v-"" = e 2 e -\t Je 



w 



s-l 



'i-py- 1 



e e 'e 

(it+a-l) 



p )s-^n) e o{^f) 



a(s-l) = -tw i(a-l)w = e -eVj&- peI &- p [0(e)+0(l)] 



Substituting the above expressions in (6.17a) and using the identities 



e ~ite~ i£ _ e -it(l-ie+0(e 2 )) _ e -it-et+0(e 2 t) 
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in order to simplify the integrand of the integral in equation (6.17a), we find the 
equation 



17TS 

e 2 



1-e- 



t hm 

£->0 

e 2 t^0 



1-f e ^P-l) e - 1E (l_ p )^lrfp 



1 _ e it(p-l} 



0<r/<t. (6.21) 



The transformation p — 1 — | maps this equation to equation (6.12). 

We note that a detailed analysis in the neighborhood of p = 1 shows that similar 
estimates are also valid near p — 1. 



A The Asymptotics of T(l — s) and %(s) 



The asymptotic formulae given in theorems 4.1 and 4.4 involve the functions x( s ) an d 



r(l — s). In this appendix, we derive the asymptotics of these functions as t — y oo. 

We use the following well known formula for the asymptotics of the Gamma func- 
tion (see for example Olver [6j, page 294): 



r(s) = e~ s e 



s „s In s 



2n 



12s 



O 



s — y oo, I arg s\ < tc — 5, 5 > 0. 



Hence, as t — > oo, 



r(s) = V2Ke- <r - it e alnm ( 1+ *) ] 12it( - 1+ ^ 



Vi 

x 



e * 



1_ Sf + 



(A.l) 



1 27TS 17T 



27rt s a e 2 e * e e 2 * 



, 1 - 6cT ^ 

1 + ^r + % 2 



1 j-KS VK_ -J. 

2irt 2 e 2 e 4 e 



1 Z7TS V7T_ -J. 

2nt 2 e <e 



^ 2 „ / 1 



1 + ^r +0 b 



< o- < 1, t — y oo, 



where c(cr) is defined by 



a 



c(a 



'I -a) 



1 

12' 



Replacing a by 1 — a in (A.l) and taking the complex conjugate of the resulting 
equation, we find 



T(l-s) 



1 -iTT(l-s) ijL .. 

2nt2 e 2 e*e 
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< a < 1, t — >■ oo. 



(A.2) 



On the other hand, the definition (2.2) of x(s) implies 



<*0Mi - .) = ^^r (s) 



e f _ e -f (1—) 



17TS 

e 2 + e 2 r(s). 



Replacing in this equation T(s) by the rhs of (A.l ) we find 

1 



X(l - s) = (2ir)2~ s t s -2e-^e 



< a < 1, t — > oo. 

(A.3) 



Equation (A.3) together with the identity (2.15) imply the following asymptotic ex- 



pression for x( s ) 



X {s) = (2n) s -k 1 " s e^e it 



< a < 1, t — y oo. 



It is of course straightforward to extend the above asymptotic formulae to higher 
order. 



B Numerical Verifications 



B.l Verification of Theorem 3.1 



Letting a = 1/2, the error term in equation (3.2) is given by 





t = 10 


t = 10 2 


t = 10 3 


7] = t 2 


-(10.4 + 5.22i) x 10" 5 


(-10.2 + 2.97i) x 10" 9 


(15.1 -4.46z) x 10- 13 


3 

7] = t 2 


-(4.00 + 4.19i) x 10~ 3 


-(1.40 + l.lli) x HT 5 


-(7.80 + 9.8M) x 10- 8 



Note that the error is proportional to t 3 /rf +a as expected. 



In order to demonstrate the effect of the higher order terms in (3.2), we also 
consider the difference between Q{s) and the first term on the rhs of (3.2), i.e. the 
difference 



cm - E 

For a = 1/2, this difference is given by 



n 



n=l 





t = 10 


t = 10 2 


t = 10 3 


7] = t 2 


-0.291 + 0.274i 


-0.341 + 0.207i 


-0.380 + 0.121* 


7] = £ 2 


0.266 + 0.0471i 


0.127 + 0.020i 


0.0360 + 0.0612i 
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Similarly, the difference between £(s) and the first two terms on the rhs of (3.2), 
i.e. the difference 



rjLl 

L 2tt J 



cm- E 



n 



. 71=1 



1 (V_ 
1 - s \2tt 



l-s 



is given by 





t = 10 


t = 10 2 


t = 10 3 


T] = t 2 


-(8.27 + 6.52i) x 10~ 2 


-(6.95 + 10.3z) x 10~ 4 


-(3.40 + 10.6i) x 10~ 4 


77 = £2 


(1.58 - 1.50i) x 10- 1 


(9.37 - 26.0z) x 10~ 3 


(-4.93 + 3.32i) x 10- 3 



The above example illustrates that the inclusion of the higher order terms in (3.2) 
improves the convergence of the asymptotic series considerably. 



B.2 Verification of Theorem 13.21 



The error term in equation (3.22) is given by 





t = 10 


t = 10 2 


t = 10 3 


(7 = 


(1.97-3.8M) x 10- 3 


(-7.76 + 65. li) x 10~ 7 


(5.62 - 3.40i) x 10" 9 


a = 1/2 


(2.23 - 4.34z) x 10~ 3 


(-2.74 + 23.5*) x 10- 7 


(6.46 - 3.82*) x 10~ 10 


a = 1 


(2.42 - 4.78z) x 10~ 3 


(-9.41 + 82.5*) x 10~ 8 


(7.13-4.22i) x 10~ n 



Note that the error is proportional to t a 3 as expected. 
B.3 Verification of Theorem 14.11 



Letting a = 1/2, the error term in equation (4.2), i.e. the term 



T(l - s)e-^r ] a - 1 x 



'0(2), 1 <T]<U< 00, 

, _ At 4 ! 

0{e ^ + |r), U < 7] < y/t < 00, 



is given by 





t = 10 2 


t = 10 4 


t = 10 6 


77 = 10 


(3.04 + 7.27i) x 10- 3 


(8.05 - 2.53z) x 10~ 6 


(84.1 - 5.12i) x 10~ 10 


rj = t 1 / 4 


(45.4 - 6.75i) x 10" 5 


(8.05 - 2.53i) x 10" 6 


(-443.6 + 6.91i) x 10" 7 


T) = t 5 / 12 


7.27- 1.45z) x 10" 1 


(-8.69 + 9.53i) x 10" 5 


(4.38 - 13.7*) x 10" 6 



In order to illustrate the increased accuracy achieved by including the higher order 



terms in (4.2), we also display the corresponding table when only the first two sums 
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on the rhs of (4.2) are included: For a = 1/2, the difference^] 
is given by 



n=l 



n=l 





t = 10 2 


t = 10 4 


t = 10 6 


7] = 10 


(5.55 - 14.1i) x 10~ 2 


(-14.4 + 7.920 x 10" 3 


(-15.9 + 4.40i) x 10" 4 




(4.76 - 7.53i) x !0" 2 


(-14.4 + 7.920 x 10" 3 


-(26.3 + 6.92i) x 10" 3 


77 = t 5 / 12 


—(2.14 + 2.15t) x 10- 1 


(3.09 - 1.910 x 10~ 2 


(-5.70 + 10.4i) x 10~ 3 



B.4 Verification of Corollary 4.3 



Letting a = 1/2 and iV = 2, we find that the error term in equation (4.22) is given by 





t = 10 


t = 10 2 


t = 10 3 


7 

7] = £l2 


-(4.08 + 2.470 x 10- 1 


(-2.39 + 3.500 x !0 _1 


(-6.42 + 15.90 x 10- 2 


.M 

7] = ti 


(6.92 + 55.70 x !0~ 2 


(9.93 + 25.70 x 10~ 2 


-(8.75 + 12.20 x 10- 3 


7] = lOt 


(1.40 - 1.120 x 10-2 


(2.24 + 5.280 x 10- 4 


(2.45 + 186+0 x 10~ 7 



B.5 Verification of Theorem 4.4 



Letting a = 1/2, the error term in equation ( |4.30 ), i.e. the term 
e- ins T(l - s)e-^7 ] a - 1 > 




1 < 7] < £3 < OO, 
+ *5 < 7] < y/t < OO, 



is given by 





t = 10 2 


t = 10 4 


t = 10 6 


n - / 27Tt 

1 y 100 


(-5.96 + 8.830 x !0" 4 


(6.44 + 1.300 x 10" 4 


(4.64 + 69.70 x 10~ 7 


77 = V2wt 


(3.59 - 10.80 x !0" 3 


(2.72 - 28.10 x 10~ 5 


(104.6 + 8.970 x 10" 7 


7] = V200nt 


(13.7-5.210 x 10~ 6 


(-9.02 + 2.270 x 10" 4 


-(14.2 + 6.740 x lO" 7 



Remark B.l In the particular case of rj = y/2irt, the formula of Siegel (Eq. (32) 
of [8]) and the corresponding formula of Titchmarsh (Theorem 4.16 of [9]) are equiv- 



alent to the formula of our theorem 4.4 This equivalence can easily be checked 



3 This difference is exactly the error in the "approximate functional equation" of Hardy and Lit- 
tlewood cf. [4]. 
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numerically — for example, for iV = 3 all three formulas yield the same numerical val- 
ues for the error terms. In this regard, we note that Titchmarsh states Siegel's formula 
in terms of the function ^(a) defined by 



„2 c , ■ ,2 



— r v e my 2 s J [e* + oiD, cos7r(^- - a - 5 

*(a) = / dw = ^ 

2tt t e w — 1 cos 7ra 



which is related to our function $(r, u) defined in (4.28) by 



#( a ) = -i$(-i }a - _j e -M^-|) ) aGC . 
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